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Iteration Procedures for the Dirichlet 
Difference Problem 


Introduction. A fundamental paper in the theory of solving partial 

differential equations by iteration has been written by FRANKEL.! In this 

paper Frankel discusses the RICHARDSON and LIEBMANN procedures and 

| their corresponding accelerated procedures. The latter are termed the 

extrapolated Liebmann and the second-order Richardson procedures. From 

his paper three disadvantages of the second-order Richardson procedure as 
compared to the extrapolated Liebmann procedure are implied. 


(1) The extrapolated Liebmann procedure appears to converge twice 
as fast. 

(2) The extrapolated Liebmann procedure seems to require fewer Opera- 
tions for each step. 

(3) The second-order Richardson procedure appears to require the 
memory of two sets of values as opposed to only one set for the extrapolated 
| Liebmann procedure. 


In the present paper all of these disadvantages will be removed for the 
case of the Laplace difference equation. The disadvantage of the simple 
Richardson procedure versus the simple Liebmann procedure will also be 
removed. Since the Richardson procedures do not require rapid access 
memory and also have simpler error eigenfunctions, the Liebmann pro- 
cedures will no longer be viewed as superior (see, e.g. SHELDON & THomas'). 

In addition it will be shown that knowledge of the optimum value of 
the parameter(s) for either accelerated procedure immediately determines 
the optimum value of the parameter(s) for the other accelerated procedure. 
These parameters will be determined if the maximum eigenvalue of the 
simple Richardson procedure is known, and some means based on this for 
estimating these parameters for regions other than rectangles will be 
presented. These regions will be bounded by vertical or horizontal lines 
connecting the equally-spaced lattice points. Curvilinear boundaries may 
be approximated by such boundaries. 

The notation of FRANKEL! will be used and some familiarity with this 
paper will be assumed. 

Some ideas presented here will be of small use if an entirely electronic 
computer is being used, as they will be expected to save only a few iterations. 
However, they will be of interest if a punched card machine is being used. 

The Richardson Method. In this method the correction process applied 
to each ¢" consists of the addition of a positive multiple of L@" (for each 
interior point). Thus 


ont = due t+ aldjn 
= 3,2 + aL ¢j—1e + O42 + Ope + Opep — 4¢52]. 
125 
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The error at each stage is 
n an n 
€j,k = $j,2 — i,k 
The error eigenfunctions of the operator L for a p X g rectangle are: 
pia F : , =1,2,---,p—1 
2 = sin (xrj/p) sin (xsk/q) {; rt iy 4 ae. 
The corresponding eigenvalues are 
L¢,2) = — 4[sin? (xr/2p) + sin? (4s/2q) ]. 


The eigenvalues Ly.) of the Laplacian operator L and the function 
t = cos (xr/p) + cos (xs/q), which is defined in the second of equations (27) 
of Frankel,! satisfy the relation, 


(A) Lo, = — 4+ 2¢. 

Therefore the eigenvalues of the iteration operator K = 1 + aL are 
(B) Kon @ 1 — 40 + 2at. 

In particular, for a = } (shown in Frankel! to be the optimum a), 
(B’) Kos) = t/2. 


For convenience denote the eigenvalues of the simple Richardson operator 
with a = }, for any region, by ¢/2. It will be seen that these eigenvalues are 
intimately connected with the eigenvalues of the other iteration processes 
and since they are the simplest they will be referred to as the “reference” 
eigenvalues. 

If ¢;,, is an error eigenfunction of the simple Richardson process (formula 
(5) of Frankel"), with eigenvalue t/2, then (—1)*+*e;, is also an error eigen- 
function with eigenvalue —#/2. Thus the non-zero reference eigenvalues 
occur in pairs with the same absolute value but opposite sign. 

It follows that if Zo and L,, are the smallest and largest in magnitude of 
the Ly,.), then 

Hue =1+ 3Lo —. [1 + 41Ln | _ (- 4tmex) 

or 


(C) Lo+ L. = — 8. 


See also MILNE, section 85. 

In Frankel! it is shown that for a = } (the optimum value) the error 
eigenfunctions of lowest and highest frequency decrease at the slowest rate. 
But there is an a which removes the latter completely, still decreasing all 
components. This a is 


(D) ao = 1/(4 + 2tmax). 


Even if tmax is unknown, an a slightly larger than } is sure to strike a heavy 
blow at this highest frequency error, as well as at other high frequency 
errors. For example, for a 20 X 20 square, ap is, to five places, 0.12577. For 
every eigenfunction corresponding to a ¢ < 0, there is an a in the range 
% <a < } such that this removes completely that eigenfunction. [Of course, 
there is an a which will remove any specific eigenfunction, but unfortunately 
when ¢ > 0 the a is such that the higher frequency errors (always present 
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due to round off) may increase tremendously. | If the problem is first solved 
using a coarse grid and then values assigned by interpolation, the error 
components will tend to be of the higher frequency type, and may be 
reduced in this manner. (See FLANDERS & SHORTLEY’, p. 1330.) 

In this paper? the lowest frequency eigenfunction has been picturesquely 
termed a “‘pillow function.” This component of the error is the most trouble- 
some to the iteration process, and the initial values should guard against it 
as much as possible. Thus assigning the average of the boundary values for 
all initial values will in general be better than assigning all zeros. It can be 
shown that there always exists a constant value for ¢9,2 which has no pillow 
function error, and the average of the boundary values is a reasonable 
approximation to this constant. This step corresponds to a block relaxation 
of SOUTHWELL.*:7 

Further block relaxations may be used (especially if the region has lines 
of symmetry) to reduce initially other low frequency error components. 

The nearness of $tmax to unity for regions with a large number of interior 
points is an indication of the condition (see TaussKy’) of the system of 
linear equations. If the remaining error is of low frequency, it may seem 
that the difference eouations are satisfied even though the approximate 
solution is far from the true solution. A more sensitive test is the following. 
It isan analogue” of the theorem on the vanishing of the integral of the normal 
derivative of a harmonic function around a closed path, and is obtained by 
summing the entire set of difference equations for the region. The result is 
that the sum of the values at the boundary points must equal the sum of the 
values at the interior points adjacent to the boundary, with the exception 
that at corners some values are counted twice and some not at all, in an 
obvious manner. The points are merely paired at opposite ends of segments 
of grid lines, and those having no mate are omitted. This principle must 
also be satisfied for all paths interior to the region, and may be used as a 
check for any low frequency error. 

Even if one desires to use the other methods described in Frankel,' the 
usefulness of the simple Richardson method is not at an end. The optimum 
a and the optimum a and 8 in these methods will be shown to be determined 
by the reference eigenvalue of the pillow function. This eigenvalue is $tmax- 

FLANDERS & SHORTLEY? give a method for determining this eigenvalue 
experimentally. One assumes zero boundary values and approximates a 
pillow function. By iteration or by the non-iteration method of this paper 
one attempts to remove the other components present and then the rate at 
which the pillow function is reducing is its eigenvalue. A good approximation 
to the eigenvalue can then be obtained by iterating once more with a = } 
and using the Rayleigh quotient formula: 








X oi ee x (es) 
tmax ~ j,k or tmax ~ j,k ‘ 
2 => GaP 2 Leen 
j,k 


A check on the results is presented here. If one actually has the pillow func- 
tion, multiplying it by (—1)*** (changing the sign at alternate points) will 
produce the eigenfunction of highest frequency, and then ap (formula (D)) 
will reduce it to zero in one iteration. 
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If one desires to solve a complete problem by the simple Richardson 
process, it may be observed that two independent computations are being 
carried on when using the optimum a of }, either of which may be omitted. 
This effectively doubles the convergence rate. The formula is 


oye = F[b3-ne + OF4ne + Oe + Oey]. 


Since ¢3,z is not used, it need not be computed. This means that ¢9,. may 
be assigned at points where j and k have, say, like parity (both even or 
both odd). Then ¢},: is computed at the points where j and & have unlike 
parity. In general, ¢3,: is computed for even at points where j and k have 
like parity, and for odd m at points where j and k have unlike parity. The 
method then becomes a Liebmann method in the sense that latest values 
are always used. 

The Liebmann Method (a = 3). In the Liebmann method the new 
value of ¢7*' is used as soon as it is computed. Thus 


oie = Oe t+ aldstie + Oipne + OF + Pep — 462]. 
With a = }, this becomes 


oye = F[(e3the + Shane + OF + Peg]. 


The Extrapolated Liebmann Method. In Frankel! it is shown that the 
optimum a (equation above) is determined by the equation (35) 


bx —4a+1=0. 


For the rectangle, tmax = cos (x/) + cos (x/q). But for other regions, tmax 
is not in general known, except that it is twice the reference eigenvalue of 
the pillow error eigenfunction. This has not been proved here, but it follows 
as a corollary of results presented in YouNnG.’ It is shown there that it is 
better to overestimate the optimum a than to underestimate it. Further 
discussion of the problem of obtaining the optimum a will be presented in 
a later section. 

The Second-Order Richardson Method. The iteration equation for the 
second-order Richardson method is: 


opt = O32 + aLdje + B(¢}.2 — oj2)- 


It appears that two lines of starting values are necessary. It is not clear 
what relation they must bear to one another or how this relation will affect 
later iterations. [If the coefficient of an error eigenfunction is A in ¢33' and 
@ in $3.2, its coefficient in ¢7%' will be a(1 — 4a + 8 + 2at) — BA. For 
B = 4a — 1 this reduces to 2aat + (1 — 4a)A ]. Actually this method as- 
sumes that the simple Richardson method has the sequence of values “‘going 
in the right direction” and therefore pushes them a little further. 
Formula (45) of FRANKEL’ has a misprint. It should read 


K* =1 — x(p* + *)!/v2. 


In the third line below this formula the word ‘‘advantage’’ should be ‘‘dis- 
advantage.’”’ This disadvantage, number (1) of the introduction, will now 
be removed, Also the above initial difficulty will disappear. 
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One of. the formulas defining the optimum a and £ in Frankel! is 
a(Lo + Ln) = — 2(1 + 8). 
By formula (C) this simplifies to 


(E) B = 4a — 1. 
The other formula in Frankel! is 

aly = — (1 — #})? 
or 


1+ aly + B = 261. 
By formula (A) 


Lo = — 4 — 2tmax- 
Thus 
1 — 4a — 2atmax + 4a — 1 = 2(4a — 1)! 
or 
(F) Ob mex — 4a + 1 = 0. 


But this is the same equation as equation (35) of Frankel.' Thus the optimum 
a in a given region is the same in the second-order Richardson method as 
in the extrapolated Liebmann method. 

Note that as tmax — 2, the optimum a — $. 

Substituting 4a — 1 for 8 and writing out the operator L in full, the 
iteration process becomes: 


Ope = Oe + aLOjte + Opps + 672-1 + O24 — 4652] 
+ (4a — 1)(¢}.2 — oF) 
or 


O7e = of 31,2 + O34.2 + O72-1 + Oheu] + (1 — 4a)opr. 


The formula has already simplified considerably. Disadvantage (2) of 
the introduction has disappeared. ¢3,, has dropped out. As in the simple 
Richardson method, it is unnecessary to compute it at all. That is, if one 
has values at alternate points, say where j + & is odd, and values in the next 
iteration at points where j + k& is even, then one can compute the values in 
the succeeding iteration at points where j + & is odd, and from these the 
values in the iteration after that at points where j + k is even again. The 
missing values are independent of the values computed, and vice versa. 
Thus only half the computation is necessary. This cancels the factor of two 
disadvantage over the extrapolated Liebmann method. 

This method of computing removes another disadvantage, number (3) 
of the introduction, of the second-order Richardson method for machines 
with limited memory. At no time is it necessary to remember more than one 
value at any point. The method is a Liebmann method in the sense that 
latest values are always used. In fact this method and the extrapolated 
Liebmann method are both special cases of the ‘‘successive overrelaxation 
method”’ developed by Young,’ depending only on ‘‘ordering’’ of the points. 
This reduced second-order Richardson method is actually produced by an 
“alternate diagonal” ordering of Young.* However, there is one difference, 
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which can perhaps be used to advantage. In one case there are ‘“‘missing”’ 
values. In the other case two successive iterations are collapsed into one. 
This of course changes the eigenfunctions. In the reduced second-order 
Richardson method, the eigenfunctions are still the simple Richardson 
eigenfunctions. 

Use of a Richardson procedure now permits the lowest frequency error 
components to be attacked both by initial block relaxations and by solution 
first with a coarse grid. The highest frequency error component can be 
eliminated in one iteration by the simple Richardson procedure with a=ap 
as given by formula (D) and other high frequency components can be 
effectively eliminated in a like manner. 

However, unlike the simple Richardson procedure, for the second order 
Richardson procedure all error components become equally ‘‘troublesome,”’ 
at least in the simple sense of asymptotic decay rate." Attacking them indi- 
vidually no longer gives such an obvious advantage. However, examine the 
expression 2aat + (1 — 4a)A (given above in brackets) for the coefficient 
of an error eigenfunction in terms of its coefficients in the two previous lines. 
Observe that one will probably obtain a second line from an initial approxi- 
mation by an iteration using the simple Richardson procedure with a = }. 
Then a = At/2 and the expression becomes (af? + 1 — 4a)A. Considered 
as a function of ¢ this is a parabola with minimum at ¢ = 0. This suggests 
that one will probably gain by making some effort to eliminate the lowest 
and highest frequency components and then iterating at least twice by the 
simple Richardson procedure with a = } before starting the second order 
procedure (assuming, of course, that the computer being used is not of 
the highest speed). 

The Reference Eigenvalue for the Pillow Function. Since the reference 
eigenvalue for the pillow function is $¢max and the optimum a in the faster 
converging methods is determined by the smaller root of the equation 


fax — 4a+1=0, 


this eigenvalue if of prime importance. A practical method for approximating 
its value, using the Rayleigh quotient, has already been mentioned. Un- 
fortunately this method in general gives a value which is too small, and as 
mentioned previously, it is much better to have a value which is too large. 
Since a formula is known for the rectangle, it is useful to observe that if the 
region in question is covered by a rectangle, the desired value is less than 
the reference pillow eigenvalue of the rectangle. Similarly a rectangle in- 
terior to the region will have a smaller reference pillow eigenvalue. In 
making these estimates, better results can be obtained if it is realized that 
for a fixed number M of interior points the reference pillow eigenvalue is a 
measure of the compactness of the region. E.g., an unbounded number of 
interior points can be strung out so that the maximum reference eigenvalue 
is less than one-half. But only four interior points arranged in a square will 
have a reference eigenvalue of one-half. This means that any arms jutting 
out from the region can be cut off and used to fill in any part of the covering 
rectangle without reducing the maximum reference eigenvalue of the region 
in question. This also means that, unless the region is almost circular, an 
upper bound to the reference eigenvalue is cos r/(1 + M*), being the 
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maximum reference eigenvalue of a square with M interior points, and a 
square being a rather compact region. 
Time Estimates. The tables of Frankel,’ page 74, can now be changed to: 


Method T 
Laplace-Richardson (a = 4) 4r 
Laplace-Liebmann (a = 3) 4r 
Laplace-Liebmann (optimum a) 6r 
Laplace-2nd-order Richardson 6r. 


Approximate Total Calculating Time. 
N=10 N=20 WN Large 
Laplace-Richardson (a = 3) 108 10s . 
Laplace-Liebmann (a = }) 5-10%r 810% 5.6N*r 
Laplace-Liebmann (optimum a) ‘ P 
Laplace-2nd-order Richardson i'r sodid did 
Here + is the mean time required for each arithmetic operation, T is the time 
for the evaluation of one ¢j,z and N is the common value of # and g for a 
square with (N — 1)? interior points. 

The change in the extrapolated Liebmann method (for optimum a) 
from 7r to 67 in the first table is because terms may be collected in the 
iteration formula to give: 


ie = aLdyte + Oypne + Opes + ORet] + (1 — dade 


and 1 — 4a may be computed when a is computed. 

The improved second-order Richardson method now has the same rate 
of convergence as the extrapolated Liebmann method, requires no more 
memory, and does not require rapid access memory. Thus it would seem to 
be a superior method for general use. 


James D. RILEY 
U. S. Naval Ordnance Lab. 
Silver Spring, Md. 
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A Method for Calculating Inverse 
Trigonometric Functions 


Introduction. The problem of calculating an angle from a known trigo- 
nometric function of that angle occurs frequently in the numerical analysis of 
physical systems. The purpose of this paper is to introduce a new method 
which is particularly suitable for the I. B. M. Card Programmed Calculator. 

Methods. There are several methods of determining an angle from its 
trigonometric functions; the more commonly employed ones are: 


1. Table look up. 
2. Infinite processes 

(a) Power series 

(b) Continued fractions.! 
3. Approximations.” 


Which of these methods is the best to use depends primarily on the com- 
puting equipment that is to be used. For desk calculating, tables are usually 
best ; for computers with limited storage capacity, the infinite processes are 
usually used ; on large scale computers, the approximations are most efficient. 
Since the IBM-CPC has a very limited amount of high speed storage, the 
common practice has been to use a Maclaurin’s series with a variable number 
of terms.* Other methods have been devised‘ that speed up the process by 
first transforming the argument into a smaller number so that the series 
will converge more rapidly. 

The New Method. The above discussion helps to illustrate the problem 
facing the programmer when programming a method for determining an 
angle from its trigonometric function. An entirely different method for 
finding inverse functions will now be described. This method is more of a 
logical process than a mathematical one. The following notation will be used: 


y = trig (A, B) 


which reads, ‘‘y is the inverse trigonometric function of A and B.” 
Furthermore, 
A 
B 


sin y, 
Cos y. 


In the case both the sine and cosine of y may be found independently, 
this method will determine the angle in its proper quadrant. If the cosine 
must be found by an identity from the sine, then only the principal value 
will be computed. 

Let us assume that the sine and cosine of some angle, y, is known and 
it is subject to the following restrictions: 


o<¥% <*¢, 
(m = 0,1, 2, «++, 25). 


mp4 
Lek 
The excluded values are special cases and will be discussed later. We first 


calculate sin 2y and cos 2y by suitable trigonometric formulae. From the 
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first restriction, 0 < y < m, we have sin y > 0; if 


sin2y>0, then 0<2y<a and O<y< 


ml a 


but if 
Ld 


sin2y <0, then w+ <2y <2 and 9 


<y<e. 


It is obvious then that a balance test of sin 2y will determine the quadrant 
of y; this information may be retained by storing the initial value of the 
quadrant. Calling the stored quantity 2, it is seen that if 


0<9¥<5, then 2, = 0, 
but if 
bal 


< hen 2, =~ 
2 y<m, then 1=% 


From sin 2y and cos 2y we compute sin 4y and cos 4y and analyze in the 
same manner. 


If sin2y>0 and sin4y>0, then 0<y<- and 22=0, 


4 
° e o Tr rT T 
if sin2y>0O and sin4y <0, then 4 <y< 3 and 22= re 
Gir « : 3 
if sin2y<0 and sin4y>0, then <y< r* and 22= z. 
: Y : 3a 3x 
if sin2y <0 and sin4y <0, then ia <y<nw and = 7° 


where <2 is the initial value of the octant in which y lies. A comparison of 
the 2»’s and 2,’s will show that when 


sin4dy>0, then 22 = 2 
and when 


sin4y <0, then 22, = 2,+ 7 


In fact, if the same analysis is extended the following generalization will 
be apparent: 
When 
sin2*y>0, then 2, = 2a-1, 
and when 


sin2"y <0, then 2, = 2ei+ “ 


The =, will be the initial value of the interval, of width #/2*, in which y 
lies. If this process is continued until 2/2" is less than the maximum allow- 
able error, then 2, will be the value of y. If the error is to be less than 
5 X 10-* then » must take on all values up to 26. 
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The cases where y = 2/2 cause difficulty because 


sin24y=0 if n2m. 


Actually, the number =, is the answer, y, when m = m since 2, = Im = =: 
Therefore, the first zero should be treated as a ‘“‘minus’”’ and all subsequent 
zeroes as “‘plus’’. 


If 
y=0 or y=rf, 
then 
sin y = 0 
and 
cosy=1 or cosy= —1; 


therefore, if sin y is zero, then a separate test must be made on the cosine 
to determine if y = 0 or y = z. 

The discussion thus far has assumed that y is in the first two quadrants; 
the method can be extended to third and fourth quadrant angles by three 
methods: 


1. The absolute value of sin y may be used and the final answer multi- 
plied by minus one if sin y was originally negative. This gives answers 
in the range — r <y Sr. 

2. If sin y is negative, set Zo = 2; in the previous discussion it was 
implied that 2» = 0. This will give answers in the range 0 S y < 2z. 

3. If sin y is negative, set 2» = — w. This puts the same limits on y 
as method 1. 


In most cases, either method 1 or method 3 is preferable to method 2. 

Method 1 is illustrated in the flow diagram, Figure 1. In practice it is 
more convenient to store c, = 2,/m and p, = 1/2” instead of 2, and 2/2", 
respectively. Since p, = $fn-1, any rounding error will be positive and the 
sum g, will include the errors of all ,’s where k < n. This error may be as 
large as six in the last digit. This error may be reduced by always rounding 
to an even number or by carrying an extra digit in », and o,; obviously 
both methods have disadvantages. The following is still another way in 
which the rounding error may be reduced and also save storage space. Let 
ad, = Oora, = 1, depending on whether the term is to be omitted or added; 
also let m = 0,1, 2, ---, N; then 


y= =n = TON 
and 


: 1 
factoring out QW? then 


1 
on = gna" + 22%? + a32"-* + --- + ay) 


or 


2% ow = {{[(2a1 + @2)2 + a3 ]2 + a4}2 +---+ ay}. 
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The quantity Zy = 2” cy may be computed by the following formula: 
Zatti aa 2Zn + Qn+1 


where Zo = 0 and m has the same meaning as before. The quantity thus 
computed is a whole number and is therefore exact, the only rounding error 
will be in the final multiplication by +/2*. The number of times Z has been 
computed is subtracted from N so that the process may be stopped when 
this quantity, ,, is equal to zero. This also saves storage space since only 
two digits of N need to be stored rather than seven or more of p,. The flow 
diagram illustrating this method is shown in Figure 2. 

In computing sin 2"*'y and cos 2"*'y both sin 2"y and cos 2”y must be 
used in each ; otherwise the condition 


(sin 2"+1y)? + (cos 2"+y)? = 1 


will not always be satisfied. As an example, suppose the following formula 
is used : 
cos 2*ty = 2(cos 2"y)? — 1. 


If sin y < .0003162, then cos y = 1.0000000 (the error being less than 
5 X 10-8) and cos 2*y = 1.0000000 for all 1; this is obviously incorrect. 
A better formula to use is 


cos 2"t1y = (cos 2y)? — (sin 2*y)?. 
By defining the quantities 


F(n, y) = cos 2"y + sin 2*y 
and 

G(n, y) = cos 2"y — sin 2*y, 
then 

cos 2"+1y = F(n, y)-G(n, y) 
and 


sin 2"t1y = F(n, y)? — 1. 


Although the above formulations are clumsy, they are necessary in order 
that the calculation can be made on the IBM 605 calculator. The pro- 
gramming for an unmodified 605 is given at the end of the article. All three- 
digit storage units are put on an 8-6 assignment and one full sweep of 
program steps is impulsed. The quantities sin y and cos y are read into 
Factor Storage 1, 2, 3 and 4 and the answer is read out of the third through 
tenth positions of the counter. During computation, F1, 2 is used to store 
sin 2"t1y and F(n, y); F3, 4 is used to store cos 2**'y; Z, is stored in G1, 2 
and 5, in G3. The following modifications would save program steps but 
would require additional tabulator wiring. General storage 1 and 2 can be 
cleared by the tabulator instead of on program step 10 to set Z) = 0. The 
quantity b) = N may be read from a card or emitted from the tabulator 
into G3. If the tube located at panel 1-7T is removed, the absolute value 
of a number may be read in on one program step instead of taking three. If 
the final multiplication is carried out on a separate card cycle, steps 50 
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through 60 may also be eliminated. This programming allows y to be any 
value in the region — x < y < x, but excludes the possibility y = x. 

In the worst cases, this method requires a little over two seconds which 
in CPC operation is less than six card cycles. This compares favorably with 
the time required by other methods to attain the same accuracy, and unlike 
most other methods does not require extra card cycles for determining the 
quadrant of the answer. 
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MORI 
Emit 6 
ECRI— 
P. U. Rpt. 
Del. 
ECRI+ 
F1i2RI 
Zero Test 


605 PROGRAMMING 


Suppress 
on 
Ri 2nd 1 
Prog. Exp No. 1 1 
Bal. Test 1 
Reset Re _ 
1, — 
Reset 
Bal. Test 
1, Non zero 
P. U. Gp. Sup. 1 
ok 
Ro 6th 
Ro 6th 
Ri 2nd 
RO 3rd 
-,2 
Non zero 
Ro 6th 
Ro 6th 
Ri 6th 
Ri 6th 
F1i2RI 
Ri 2nd 
Ri 5th 
Ro 6th 
Ri 6th 
Ri 6th 
Ro 3rd 
Ri 6th 
Prog. Rpt. Zero 
Ri 2nd Non zero 
Ri 2nd Non zero 
Prog. Exp 3 Non zero 
Ri 5th Non zero 
Ro 2nd Non zero 
Non zero 
Emit 6 Non zero 
Ro 6th Non zero 
Non zero 
Non zero 
Ri 2nd Non zero 


Notes 


Program expansion No. 1 
is G3RI, Ri 6th 


* D.O. Rpt. Del. if Repeat 
selector one is up. 


Program expansion 2 is 
ECRI-+, Ri 6th 

Program expansion 3 is 
Ri 3rd, Ri 4th 


* Multiply plus if repeat 
selector 1 is normal, 
minus if transferred 


* Multiply plus if repeat 
selector 1 is normal, 
minus if transferred 
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Steps Description of Operation 
8, 11-13 Compute cos x — sin x and store 
14-16 Compute cos x + sin x and store 
17-23 Compute cos 2x and store 
24-27 Compute Z,,,; and store 
29-45 Compute sin 2x and store 
46-48 Modify 6, 
50-60 Multiply Zy by 7/2” 


RoBerT L. LAFARA 
U. S. Naval Ordnance Plant 


Indianapolis 18, Indiana 


on TEICHROEW, “Use of continued fractions in high speed computing,”” MTAC, v. 6, 
p. 127-133. 


2 C. HasTINGs, Approximations in Numerical Analysis. (Form 15’s). Rand Corporation, 
Santa Monica. Sheets 8-12, 35-37, 39. 


’W. P. Hersinc, “An eight-digit general-purpose control panel,’’ IBM Tech. News 
Letter No. 3, 1951. 


4D. W. Sweeney, “A Model II four address floating-decimal coding system,”” IBM 
Tech. News Letter No. 5, 1953. 


A Partitioning Method of Inverting Symmetric 
Definite Matrices on a Card-Programmed 
Calculator 

Although certain standard lines of approach seem to exist in all methods 
of inverting matrices, the problem of inversion has not yet been solved with 
finality, and modifications on ‘‘well-worn’”’ solutions often prove to be of 
value. This paper discusses a partitioning method for inverting matrices 
which is particularly suited to the International Business Machines’ Card 


Programmed Calculator, (C.P.C.), and subsidiary equipment. 
The general problem is to invert an m“ order matrix 


_f{a _b 
u=(¢ 3) 
where a and d are each square submatrices of order m; and n2[m, + nz = nj. 
Let the inverse of M be 
“a? oe 
w= (5 5) 


partitioned in the same manner as the original matrix. Post-multiplying @ 
by its inverse, we have 


aA + 6B $13) -(] +4 
cA+dB cC+dD/ \0 1 


or in equation form 


aA +b6B=1 (identity matrix of order m;) 
cA +dB=0 (zero matrix m2 X 1) 
aC+bD=0 (zero matrix m X m2) 
cC+dD=1 (identity matrix of order m2). 








140 A METHOD OF INVERTING SYMMETRIC DEFINITE MATRICES 


Solutions of these equations are 


D = (d — ca)“ 
C= —a"bD 
(1) B= — Dea" 


A =a — aB. 


Thus to find the inverse M-', the procedure is to determine a~! from 
which we can derive D, C, B, and finally A. If the original matrix is sym- 
metric, equations (1) simplify to 


D = (d — b’a™b)“ 
(2) B’=C=-—-a"bD 
A a — adB. 


The restriction to definite matrices M assures us that (d — b’a—) = 0. 

Let b and ¢ be a column and row vector, respectively, and d a scalar; 
then the solution of (2) is obtained relatively simply once a— is known. (Note 
that D is merely the reciprocal of a scalar.) a~! can be found by partitioning 
a in an analogous manner and a™, which is of the order » — 1, depends on 
the inverse of an » — 2 submatrix. The process carried to its completion 
implies that the inverse of the original matrix M is a function of the inverses 
of smaller and smaller submatrices. Thus a well known method for inverting 
matrices is suggested: start by inverting a small submatrix, say of order 2; 
from the 2 X 2 inverse, derive the 3 X 3 inverse, and continue until the 
desired » X m inverse is computed. This method turns out to be efficient and 
very simple to handle on an automatic computer. 

The programming suggested assumes the following C.P.C. set-up: (a) 
Products can be accumulated in the 605 Electronic Calculating Unit 
counter; e.g., if a = x.y; is derived on card p and B = xey-2 is derived on 
card p + 1, then on card p + 1, a + 8 can be accumulated in the counter. 
In other words, the C.P.C. must be able to perform both a multiplication 
and an addition on one card.! (b) A second required feature is that a card 
may be read from either a normal or an alternate instruction field. 

The procedure employs (2m — 1) storage locations. 

The manual operating scheme is: 


1. Load the 527 Summary Punch with cards pre-punched in their 
instruction fields. There are m stages and at the r** stage, the 527 summary 
punches r? cards which are fed back into the 418 Tabulator at stage (r + 1); 
the r? cards are the r** order inverse a~. 

2. At the beginning of the r** stage, the operator removes (r — 1)? cards 
summary punched on the (r — 1)** stage and runs them through the collator 
which assembles them in row-by-row sequence. These (r — 1)? cards are 
manually put in the program deck for stage r. The first part of this deck is 
run through the tabulator; when the (r — 1)? cards have gone through the 
tabulator once, they are removed and placed at the end of the program deck 
for the r** stage and the deck finishes its run. A new set of r? cards will 
have been summary punched and the operator is ready to go on to stage 


(r + 1). 
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The table below gives the functional computing procedure at the r™ 
stage. The following notation is used: x is a storage bank having n locations 
[x+1,---,x-+m] and y a storage bank having m» — 1 locations, 


Ly + 1,---,9 + (# — 1)]. 


Table 1. Symmetric Matrix—Stage r 


Put into Number 
Storage of 
Step Function Location Operation Cards 
1 Load bd x+1 (Add) r—1 
x+ (r —1) 
Load d x+r (Add) 1 
2 Form ab y+1 Row by row, the (r — 1)* cards contain- (r — 1)? 
: ing a are put through the tab. The 
. first instruction field is punched so as 
y + (r — 1) to give the product of each element of 
a row of a~, from the card, times the 
corresponding element of b, from stor- 
age x, accumulate these products for 
each row and store the (r — 1) sums 
in y. 
3 Form Keep in 605 Call out the element of b from storagex, r— 1 
(d — b’a~b) Counter negative multiply by the elements of 
a~'b in storage y, and accumulate the 
products with d which was called in on 
the first card from x + r. 
Take Division 
1 
G-veu 7? *t" 
4 Calculate Summary Call out each a~'d from storage y, nega- r—1 
-—-a)D=C Punch tive multiply by D and punch out. 
5 Calculate Summary Same as 4. r—1 
—a'bD=B Punch and 
Store: 
x+1 
x + (r — 1) 
6 D Summary Take D from storage x + r and summary i 
Punch punch. 
7 Calculate Summary Run the (7 — 1)? cards through the tab (r — 1)* 
a*—a¥B=A Punch again, this time using second instruc- 


tion field. As each row of a goes 
through, take the element of a~'d from 
storage y, and negative multiply it by 
the element of B in storage x, add it to 
the element of a“ from the card- 
summary punch the answer. 


Total number of cards on the r stage. 2r°+ 1 


Notice the order of punching out: we first obtain the last element in each 
row of a~ for the first (r — 1) rows (C); then we obtain the entire r row 
(B, D); and finally the first (r — 1) rows less their last element (A). If we 
break this deck into a deck of C, B, and D (already in order) and a deck of A 
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Table 2. Approximate C.P.C. Running Time for Symmetric Matrices 


Order | oP 
n 10 16 20 25 30 35 40 pu 
number of 816 3112 5928 11,344 19,332 30,392 45,024 
cards be 
running 14 52 1 hr. 3 hrs. 5 hrs. 8 hrs. 12 hrs. me 
time min. min. 35 min. 10 min. 20 min. 20 min. 30 min. al 
we can put one deck into each feed of the collator and merge them into in 
proper row-by-row order. 

The approximate total number of cards? passing through the C.P.C. is Th 
3n* + $n’; the average speed of the C.P.C. for this method is 60 cards/minute. Sa 

If a high speed automatic calculator is used where storage is less limited, on 
then the major time consuming factor, summary punching, can be eliminated o 
and the relative computing time shortened. 

In addition to the speed of the method, the procedure has several other lat 
recommending features: if inverses of submatrices are needed, these are P- 
available by-products of the calculations. Similarly the programming of an @ 
n X n inverse includes the programming of the inverses of smaller matrices. 

An example of a 5 X 5 matrix inverted by this procedure is = 

2. «2. «| = 6 —4 1 0 0 Ri 
2 5 8 fi1 14 —-4 6-4 1 O pa 

M=1|3 8 14 20 26) Mo= 1 -—4 6 —4 1 
4 11 20 30 40 0 1 —4 5 2 at 
hy 5 14 26 40 55 0 0 1 2 1 a 
: The procedure gave the exact inverse in less than .75 of a minute. nN 


Under the assumption that in general it takes m* operations to invert an 
n* order matrix, H. HoTELLING? asserts, ‘In partitioning a matrix [for D 
" inversion ], there is an advantage in dividing the rows into equally numerous Y 
groups, since when the sum of two positive numbers is fixed, the sum of 
their cubes is a minimum when they are equal.” If Hotelling’s inference 
were correct, the method presented in this article of partitioning by single 
rows and columns would be inefficient. But in actuality the number of calcu- 
lations required to invert any matrix by partitioning is invariant under the 
mode of partitioning: 





With respect to computing time, it is most efficient to use single row and te 
column partitioning since this method enables us to store 6 and d, and also p 
at later steps a—'b, D, and B. If b were not a vector and d not a scalar, then n 
storage of these quantities would become impractical and the speed of t! 

Table 3 3 

Function Approximate Number of Operations 

a ne n 
ab nynz | t 
ca nine N 
cab ¥ nin? c 
(d — cad) ns Vv 
— a bD nn? t 
— Dea“ nn 
a —a"bB n\n | 





Total ny + 3nin2 + 3nin? + no = (m1 + m2)? = n® 


w- CO be 


— 
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operation would be slowed down further by more intermediate summary 
punch operations. 

It can be shown that if mo is the order of the largest matrix which can 
be inverted given the storage capacity, and m» < nm, then the most efficient 
method of inverting by partitioning requires inverting of, say, k submatrices, 
where the first ( — 1) matrices inverted are of order mo, and the k® 
inverted matrix is of order n — (Rk — 1)mo. 

Harvey M. WAGNER 
The RanpD Corporation 
Santa Monica, Calif. 
and 
Stanford University 
Stanford, Calif. 

This article is contained in H. M. WAGNER, “‘Matrix Inversion on an Automatic Calcu- 
lator by Row and Column Partitioning,’’ P-417, The RAND Corporation, July 14, 1953; 
P-417 also discusses in detail (1) the inversion by partitioning of non-symmetric matrices; 


(2) optimal partitioning; (3) inversion of m** order matrices where m exceeds storage 
limitations. 

1W. OrcHarpD-Hays, “The Dupiex System for IBM’s Model II CPC, A Fast Four 
Address, Double Operation, Floating-Decimal Set-Up,”” RM-1044, The RAND Corporation, 
February 23, 1953; H. M. WaGner, “‘Coder’s Manual for Duplex System” and “Stanford's 
Revised Duplex System—Wiring Manual,” Technical Reports 1 and 4, Stanford Com- 
putation Center. These references give set-ups which meet the two operations per card 
requirement. 

2 The formula and table following do not take into account any saving in cards and time 
at the beginning of the procedure when r is small enough to invert without any summary 
punching. E.g., with 19 storage locations, a 5 X 5 matrix can be inverted without any 
summary punching in 90 cards at 150 cards/minute. In general, if S is the number of 
storage locations available, then all m** order symmetric matrices for which » < N where 
jae RM ~ Le S can be inverted without any summary punching. 


3H. HorTe..inG, “Some computational devices,” Chapter X in Statistical Inference in 
Dynamic Economic Models, Cowles Commission Monograph 10, T. C. Koopmans, ed. New 
York, 1950, p. 323-328. 


Polynomial Approximations to Elementary 
Functions 


The increasing use of high-speed computing machines has revived in- 
terest in the approximation of functions of a real variable, particularly by 
polynomials. An orthodox table of function values at equidistant arguments 
may require considerable storage space in an electronic machine. In contrast, 
the coefficients of a polynomial which represents the function to a desired 
accuracy over a specified range may require very little storage space, and 
simple instructions will suffice for the evaluation of the polynomial. 

If a function is bounded and continuous in a given finite range of argu- 
ment, a powerful polynomial approximation is usually obtained by trunca- 
tion of the infinite expansion of the function in Chebyshev polynomials. 
Many properties of these polynomials, and the means by which the coeffi- 
cients in the infinite expansions can be derived, are given by LaNczos,' 
whose notation for the Chebyshev polynomials in the range 0 < x < 1 will 
be used here 

T,* (x) = cosn@, cos@ = 2x — 1. 


In this note, tables are given which ease the calculation of polynomial 
approximations to some common functions. Coefficients in the infinite 
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Chebyshev expansions valid in the specified ranges are given to nine decimal 
places, and formulae are also given for the Chebyshev polynomials. 

A polynomial approximation to any of the functions considered is ob- 
tained by truncation of the infinite Chebyshev series. The retained coeffi- 
cients are then rounded-off, and the polynomial rearranged in powers of the 
independent variable. 

As an example, an approximation to sin $x in the range — 1 < x < Lis 


5 
x mm A,T,* (x), 
n=0 
where the coefficients are given in the tables. Substituting their values gives 


sin 4x = x{1.27627 8962 — 0.28526 1569(2x? — 1) 
+ 0.00911 8016(8x* — 8x? + 1) 
— 0.00013 6587(32x* — 48x* + 18x? — 1) 
+ 0.00000 1185(128x* — 256x* + 160x* — 32x? + 1) 
— 0.00000 0007(512x" — 1280x® + 1120x* — 400x* + 50x? — 1)} 
= 1.57079 6326x« — 0.64596 4102. 
+ 0.07969 2704x* — 0.00468 1984x7 
+ 0.00016 0640x® — 0.00000 3584x"" 


In this case, the rapidity of convergence of the Chebyshev series indicates 
that the truncation error will be small, and most of the error in the above 
approximation arises from the rounding-off of the coefficients A,. Thus the 
maximum error cannot exceed three units in the ninth decimal place. 

If an approximation is required to less than nine decimal accuracy, 
unwanted rounding errors can be avoided by retaining one or two extra 
decimals in the coefficients. 

For instance, suppose that an approximation to In (1 + x) is required 
in the range 0 < x < 1, accurate to three decimal places. Examination of 
the coefficients A, in the series 


In (1+ x) = > A,T,* (x) 
n=0 
given in the tables, shows that |A,| < 0.0005 when m > 3, and that 


> |A,| = 0.000503: - -. 


n=4 


Hence 


In (1 + x) = 0.37645 + 0.34315(2x — 1) — 0.02944(8x? — 8x + 1) 
+ 0.00337(32x? — 48x2 + 18x — 1) 
= 0.00049 + 0.98248x — 0.39728x? + 0.10784x3 


with a maximum error not exceeding 0.00053. 
C. W. CLENSHAW 
National Physical Laboratory 
Teddington, Middlesex 
England 


This note is published with the permission of the Director of the National Physical 
Laboratory. 
1 Tables of the Chebyshev Polynomials S,,(x) and C,(x). NBS Applied Math. series 9, 1952. 
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Tables of Chebyshev coefficients 


sin 44x cos $4x 
nN A. n A, 
0 1.27627 8962 0 0.47200 1216 
1 —0.28526 1569 1 —0.49940 3258 
2 0.00911 8016 2 0.02799 2080 
3 —0.00013 6587 3 —0.00059 6695 
4 0.00000 1185 4 0.00000 6704 
5 —0.00000 0007 5 —0.00000 0047 
sin ax = x > A,T,*(x’), cos ax = > A,T,*(x*), 
n=0 n=0 
—-1<s$2 <1. —is2zS1. 
arctan x arcsin x, arccos x 
n ro” nN A, 
0 0.88137 3587 0 1.05123 1959 
1 —0.10589 2925 1 0.05494 6487 
2 0.01113 5843 2 0.00408 0631 
3 —0.00138 1195 a 0.00040 7890 
4 0.00018 5743 aa 0.00004 6985 
5 —0.00002 6215 5 0.00000 5881 
6 0.00000 3821 6 0.00000 0777 
7 —0.00000 0570 7 0.00000 0107 
8 0.00000 0086 8 0.00000 0015 
9 —0.00000 0013 9 0.00000 0002 
10 0.00000 0002 
arctan x = x >> A,T,*(x’), arcsinx = x >> A,7,*(2x’), 
n=0 n=0 
—-i<x<1. — 42 < x < 4v2. 
for |x| > 1, use arccos x = 44 — x >> A,T,*(2x*), 
arctan x = 4” — arctan (1/x) —_ 
0<x < v2. 


For $v2 < x < 1, use 
arcsin x = arccos (1 — x?)!, 
arccos x = arcsin (1 — x*)!. 


e* e* 
n An n A. 
0 1.75338 7654 0 0.64503 5270 
1 0.85039 1654 1 —0.31284 1606 
2 0.10520 8694 2 0.03870 4116 
3 0.00872 2105 3 —0.00320 8683 
4 0.00054 3437 4 0.00019 9919 
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5 0.00002 7115 5 —0.00000 9975 
6 0.00000 1128 6 0.00000 0415 
7 0.00000 0040 7 —0.00000 0015 
8 0.00000 0001 


eo-=m > AsT.*(z); 0< x <1. o*= >, A,7,*(s%), O< x <1. 
n=0 n=0 





log (1 + x) r(1 + x) 

n me n A, 
0 0.37645 2813 0 0.94178 5598 
1 0.34314 5750 1 0.00441 5381 
2 —0.02943 7252 2 0.05685 0437 
3 0.00336 7089 3 — 0.00421 9835 
4 —0.00043 3276 4 0.00132 6808 
5 0.00005 9471 5 —0.00018 9303 
6 —0.00000 8503 6 0.00003 6069 
7 0.00000 1250 7 —0.00000 6057 
8 —0.00000 0188 8 0.00000 1056 
9 0.00000 0029 9 —0.00000 0181 
10 —0.00000 0004 10 0.00000 0031 
11 0.00000 0001 11 —0.00000 0005 
12 0.00000 0001 

log (1 + x) =  AnT,*(x), T(1 +2) = D AnT,*(x), 
n=0 0=n 
0s 2< 1. e442 < 1. 
Jo(x) Ji (x) 

1 An n A, 
0 0.03154 0613 0 0.06942 43523 
1 —0.21461 6183 1 —0.11557 79057 
2 0.00433 6620 2 0.12167 94099 
3 —0.26620 3654 3 —0.11488 40465 
4 0.30612 5520 4 0.05779 05331 
5 — 0.13638 8770 5 —0.01692 38801 
6 0.03434 7540 6 0.00323 50252 
7 —0.00569 8082 7 — 0.00043 70609 
8 0.00067 7504 8 0.00004 40991 
9 —0.00006 0947 9 —0.00000 34583 
10 0.00000 4309 10 0.00000 02172 
11 —0.00000 0246 11 —0.00000 00112 
12 0.00000 0012 12 0.00000 00005 
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rc) x? r) x? 
Jo(x) = £47. (5), Ji(x) = 2¥ aT. (5), 
—-10 <x < 10. —-10<*< 10. 
T,* (x) 
1 
2x —1 
8x? — 8x + 1 


32x8 — 48x? + 18x — 1 
128x* — 256x* + 160x? — 32x +1 
512x5 — 1280x* + 1120x* — 400x? + 50x — 1 
2048x* — 6144x5 + 6912x* — 3584x* + 840x? — 72x + 1 
8192x7 — 28672x* + 3942425 — 26880x* + 9408x* — 1568x* + 98x — 1 
32768x* — 1 31072x7 + 2 12992x* — 1 80224x5 + 84480x* — 21504x* 
+ 2688x? — 128% + 1 
9 1 31072x® — 5 89824x* + 11 05920x7 — 11 18208x* + 6 58944x5 
— 2 28096x* + 44352x* — 4320x* + 162x — 1 
10 5 24288x" — 26 21440x® + 55 70560x* — 65 53600x’ + 46 59200x* 
— 20 50048x5 + 5 49120x* — 84480x* + 6600x* — 200x + 1 
11 20 97152x" — 115 34336x" + 273 94048x° — 367 65696x* 
+ 306 38080x’ — 164 00384x* + 56 37632x5 — 12 08064x* 
+ 1 51008x* — 9680x* + 242x — 1 
12 83 88608x" — 503 31648x" + 1321 20576x" — 1992 29440x*° 
+ 1905 13152x* — 1203 24096x7 + 506 92096x* — 140 57472x5 
+24 71040x4 — 2 56256x* + 13728x? — 288x + 1. 
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1202[A,P].—M. L. CLInnick, (a) Gear Ratios No. 43. (b) Gear Ratios No. 59. 
Privately printed, 3211 School Street, Oakland 2, California, 1953. Each 
book has 84 unnumbered pages 8.9 X 11.4 cm. and 8.26 X 14.0 cm. 
respectively, photo-offset from typescript. Price $1.00 each. 


These pocket-sized tables are designed for use in selecting appropriate 
sprocket gears in motorcycle racing events. They are triple entry tables 
giving 2D values of 

R = kr/(ec) 


for c = 10(1)23, e = 15(1)25, r = 46(1)75, k=43, 59. Unrealistic values 
of R > 15 are omitted. In the intended application r, e and c are the numbers 
of teeth in the rear, engine, and countershaft sprockets respectively. The 
clutch sprocket is assumed to have 43 teeth or 59 corresponding to certain 
popular British and American makes of motorcycle. Instructions are 
given in (b). 

D. H. L. 


1203[A,B,P].—J. K. Lyncu, Kilocycle-Radian Frequency Conversion Tables. 
Commonwealth of Australia, Postmaster-general’s Dept., Research 
Laboratories, Report No. 3726. Melbourne, 1953, 24 p. 20.2 XK 25.4 cm. 


This table gives 6S values of 
w = 20007 F 
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and w’, w and w~? for 
F = .1(.1)100(10)1000, 10* for k = 6(1)11. 


This allows the electrical engineer to pass easily from a frequency F in 
kilocycles per second to the corresponding angular velocity w in radians per 
second, or to its square, reciprocal or reciprocal square. 
For example the inductance L and capacitance C for resonance of 
frequency F are related by L = w?/C. 
D. H. L. 


1204[A,D].—G. E. REYNoLDs, Conversion Tables of Tangents or Cotangents 
to Sines and Cosines of Three Decimals. Air Force Cambridge Research 
Center, Technical Report 53-29, Cambridge, Mass., 1953. 26 p. 18.4 
X 27.3 cm. 


This is a short handy table of sines and cosines as functions of tangents 
or cotangents. Values of sines and cosines are given to 3D as consecutive 
multiples of 10-*. Corresponding 4S values of tangents and cotangents are 
listed alongside. The table is intended to be entered via the more rapidly 
moving functions tan and cot to read out corresponding values of sine and 
cos without interpolation. 

The table is in fact a table of r/(1 — r*)#, r — .0005 = 0(.001)1. Appli- 
cations are mentioned to the layout of parallel plate surfaces and the path 
of a milling cutter. 

D. H. L. 


1205[B,F].—D. R. Kaprexar, Cycles of Recurring Decimals, v. 11. (From 
N = 167 to 213 and many other numbers.) Khare Wada, Deolali, India, 
1953. Published by the author, iv + 47 p., 24.3 K 16.7 cm. Price 
6 rupees. 


This is an extension of v. I, reviewed in MTAC, v. 7, p. 238, to the extent 
mentioned in the title. Besides the table of cycles there are other tables 
as follows. 


P. 25. Table of 2" and 5” for m = 1(1)33 
P. 36-41. Table of exponents of 10 (mod p) for all primes p < 13709. 


This table is taken from a previous table of KRAITCHIK.' This reprint 
contains three errata. 


p. 36 For P = 669 read P = 659 

p.40 P = 9619, for C = 8 read C = 3 

p.41 For P = 12901 read P = 12907 

P. 43-47. Table of factors of (10° — 1)/p where e is the exponent of 
10 (mod ») for 23 primes p mostly less than 100. In connection with 
p = 47 the author claims the discovery of the prime 
193423597678916827853. 
However on the next page two factors of this number are given, 
namely 2531 and 549797184491917, the last being misprinted. 
D. H. L. 


1M. Kraitcuik, Recherches sur la Théorie des Nombres, v. 1. Paris, 1924,"p. 131-145. 
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1206[C,D ].—Akademiia Nauk SSSR. Institut Tochnoi mekhaniki i vychi- 
slitel’noi techniki. Matematicheskie Tablifsy. Desiatiznachnye Tabliisy 
logarifmov kompleksnykh chisel i perekhoda ot dekartavykh koordinat k 
poliarnym. Tablitsy funkisit. [Ten place tables of logarithms of complex 
numbers and of the transformation from cartesian to polar coordinates. 
Tables of functions] In x, arctg x, 4 In (1 + x’), (1 + x*)#. Moscow, 
1952. 116 p. 17 X 26 cm. 10.8 roubles. 


There are four 10D tables each at interval .001 and with A?: 


A. Inz,1i<2z < 10; 

B. tin (1+ 2),0<x<1; 
C. arctanx,0 <x <1; 
D. (1+ 2x)3,0<x<1. 


From such tablesInz=In (A +2B) =In A+$1n [1+ (B/A)*]+éarc tan (B/A), 
|B| <|A|, r =|A|[1 + (B/A)*], 6 = 4x — arctan (B/A) may be found; 
similar expressions exist for |B| >|A]|, since In (A+2B) =In i+In (B—iA). 

Illustrative numerical examples in the use of the tables are given on 
pages 4-6, and there is an errata slip with 15 corrections in the following 
tables. An interpolation sheet is in a pocket of the volume, of which there 
were 3000 copies in the edition. 

R. C. ARCHIBALD 

Brown University 
Providence, R. I. 


1207[D ].—H. E. Sazer, “Radix table for obtaining hyperbolic and inverse 
hyperbolic functions to many places,” Jn. Math. Phys., v. 32, 1953, 
p. 197-202. 


These tables are intended to give data from which 18D values of hyper- 
bolic functions and their inverses can be found with an ordinary desk calcu- 
lator. A similar table for circular functions has been given by the author 
in MTAC, v. 5, p. 9-11. As in the previous radix table, the author chooses 
the inverse tangent function. In this case it is 


Arc tanh (k-107) k= 1(1)9, A= 1(1)6. 
These values are given to 20D. He also gives 20D values of 


tanh A for A = 1(1)24 
and 
41Inx for x = 1.2(.2)2(1)10. 


By using the addition theorems for tanh and arc tanh and approximations to 
these functions for small values of the argument any desired value can be 
built up. Of course, the other hyperbolic functions and their inverses can be 
expressed in terms of the tangent and its inverse. 

D. H. L. 


1208[F ].—Ove Hemer, “Note on the diophantine equation y? — k = x,” 
Arkiv fir Matematik, v. 3, 1954, p. 67-77. 


This paper contains corrections and additions to the author's dissertation 
[MTAC, v. 7, p. 86]. 
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Tables 1 and 2 of this dissertation are reprinted here with the corrections 
called for in our previous review. Table 1 is now supposed to be complete 
in the sense that it gives all solutions of the equation y? — k = x* for 
0 < k < 100, whenever the equation is soluble, together with the number 
N of solutions. The largest N = 8 is for the well known case of k = 17 
discussed by MorDELL, Table 2 for — 100 < k < Ostill contains 22 incom- 
plete cases. For these, the coefficients and solutions of the corresponding 
cubic forms are given. Table 3 of the dissertation is not reproduced here. 
Table 4 is slightly enlarged giving the discriminant —D, the fundamental 
rings and units in the cubic field corresponding to all square free values of 
1 < k < 50, and for the values of k corresponding to soluble equations for 
50 < k < 100. There still remain 7 cases in which the unit is “‘not definitely 
proved to be fundamental.” 

EMMA LEHMER 
942 Hilldale 
Berkeley, Calif. 


1209[ F ].—D. H. LEumMer, Emma LeumMer, & H. S. VANDIVER, “An applica- 
tion of high speed computing to Fermat’s Last Theorem,” Nat. Acad. of 
Sci., Proc., v. 50, 1954, p. 25-33. 


Kummer defined a prime / to be regular if it does not divide any of the 
first (1 — 3)/2 Bernoulli numbers and showed that 


st+y=2', 1>2 


is impossible in non-zero integers if / is a regular prime. In this paper, the 
irregular primes less than 2000 are tabulated, together with the accompany- 
ing least prime p of the form 1 + &i, as well as three other associated con- 
stants; 2a, Qa, Qa*. On the basis of this table not only is Fermat’s Theorem 
established for / < 2000, but data are given which will greatly simplify and 
facilitate the study of units and ideals in cyclotomic fields. 

B. W. JoNEs 
Univ. of Colorado 
Boulder, Colo. 


1210[G,K ].—F. N. Davin & M. G. KENDALL, ‘“‘Tables of symmetric func- 
tions, Part IV,”’ Biometrika, v. 40, 1953, p. 427-446. 


In this set of tables the authors continue a series of fundamental sym- 
metric function tables for functions of weight not exceeding twelve. When 
the series is complete the user will be able to relate any two of the following 
kinds of symmetric functions of x;: 


(U) unitary, or elementary, a, = >> x1%2-+ +x, 
(S) sums of like powers, S, = >> x1" 

(M) monomial, such as (3, 2,1,1) = }> x:3xo2xgx,4 
(H) homogeneous product-sums h, generated by 


IT (1 — xa) = & he 


The first three parts! have already been published and are described in 
RMT 769 and 1020 [MTAC, v. 4, p. 146, v. 6, p. 224-5, see also corri- 
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gendum, v. 8, p. 188]. These relate S and M, M and U, and U and H 
respectively. The present Part IV relates M and H. Thus one may use the 
table of weight 3 to read such results as 


hyhe = (3) + 2(2, 1) + 3(1, 1, 1) 
and 
(2,1) = — 2h? + Shyhe — 3hs. 


In the arrangement of the tables advantage is taken of Hirsch’s law of 
symmetry which states that the coefficient of (p, g, 7, ---) in the expansion 
of hah.--- is the same as the coefficient of (a, b, c, ---) in the expansion 
of h,h,h,- --. Thus although each expanded symmetric function of weight w 
requires w terms each coefficient may be used twice. The whole set of tables 
for Part IV occupies no more space than the previously published parts 
which are triangular in nature. The printing is in the same small but elegant 
type of the earlier tables. 

Part V which will relate symmetric functions of types U and S will 
complete the set since a Part VI relating H and S would be identical with 
Part V except for obvious sign changes. 

D. H. L. 


1F. N. Davip & M. G. KENDALL, “‘Tables of symmetric functions Part I; Parts II and 
III,”’ Biometrika, v. 36, 1949, p. 431-449; v. 38, 1951, p. 435-462. 


1211[K].—R. S. Burtncton & D. C. May, Jr., Handbook of Probability 
and Statistics with Tables. Handbook Publishers, Sandusky, Ohio, 1953. 
ix + 332 p., 14.3 X 20.3 cm. Price $4.50. 


This book is principally a rather complete and up-to-date handbook of 
statistical methods and theory which will prove very useful; reviews con- 
cerned with this aspect will be found elsewhere. Scattered through the text 
are some 17 tables and at the end 24 tables and a nomogram occupy 62 pages. 

The tables in the text, mostly short, begin with four giving ordinates, 
areas, and percentage points for the normal frequency function. The follow- 
ing four are devoted to the normal bivariate frequency function. The first, 
if this function be written f(x, y) = k exp (—c*/2), gives c to 4S for the P% 
probability ellipse for P = 25, 50, 75, 90, 95, 99. The next three deal with 
the circular case, giving radii co of P% probability circles with c to 4S for 
P = 25, 39.3, 50, 54.4, 75, 90, 95, 99; the radii R; to 4D in standard units 
of a circular disk centered at the means over which the probability integral 
is P for P = 0(.05)1 as well as the integrand evaluated at R, to 4D; and 
values of the same probability integral to 3 or 4D over a circular disk of ra- 
dius R/o centered at the distance d/o from the means for R/o = .1(.1)1(.2)3 
and d/o = 0(.1)3(.2)6. The next two tables are concerned with the trivariate 
normal frequency function with zero correlations, the first being the exact 
analog of the first bivariate table and the second, for the spherical case, the 
radius co (c to 3D) of the sphere within which the total frequency is for 
pb = .25, .5, .75, .9, .95, .99. 

All but one of the remaining tables in the text are related to sampling 
distributions. Those for which credit is given to other sources are here 
designated by (R). If s is the standard deviation in a sample of m from a 
normal universe with standard deviation o, the next table (R) gives 6 
for which P(s > bc) = P to 3D for nm = 5(1)30 and P = .01, .05, .95, .99. 
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If r = r,Vn — 1 where rz, is the coefficient of correlation in a sample of m) 
from an, uncorrelated bivariate normal universe, the following table (R 
gives values of r; to 3D for which P(r = r;) = P for P = .001, .01, .05, .1 
and m = n — 2 = 1(1)30(5)50(10)100, 120, ©. Then comes a table (R) 
concerned with the distribution of ranges, R, in samples of m from the same 
universe. Values of E (2) are given to 4S for m = 2(1)20, to 3S for m = 30, 
50, 75, 100, and to 2S for m = 150, 200 and op,, is given to 3S for m = 2(1)20. 
Also R’/e is given to 2D for which P(R < R’) = P for P = .001, .005, .01, 
025, .05, .1, .9, .95, .975, .995, .999 and m = 2(1)20. Next the probability 
p that the fraction 4 of any continuous universe is included in the range of 
a sample of » drawn from it is tabulated to 2S (more in a few cases) for 
h = .8, .9, .95, .99 and m = 5(5)30, 40, 50, 75, 100. The next pair of tables 
(R) give equal-tail confidence intervals for the probability of success p in 
a binomial universe from which x trials have yielded S successes. The end- 
points are given to 2D for m = 10, 15, 20, 30, 50, 100, S = 0(1)50 and the 
confidence coefficients .95 and .99. Then follows a table of Stirling’s numbers 
of the first kind, S;", for m = n(1)10 and 7 = 0(1)9. The final table (R) in 
the text gives the numerical coefficients for estimating center lines and 
control limits for X and o¢ or R charts in quality control, both for process 
means and standard deviations known or unknown. These are given to 3 
or 4S for samples of m = 2(1)10(2)20 for averages and standard deviations 
and m = 2(1)10(2)14 for ranges. 
The section devoted to tables (p. 247-309) include: 


Table I. (*) p*(1 — p)** to 4D for p = .05(.05).5 and m = 1(1)20. 


Table II. > (*) p?(1 — p)** to 4D for the same values of p and n 


zrZ=zr 


as in Table I. 

Table III. » for which the incomplete 8-function, J,(x,n — x +1) = .005, 
01, .025, .05, .1 for x = 1(1)15, 20, 30, 60, © and m — x + 1 = 1(1)6, 10, 
15, 20, 30, 60. 


Table IV. (*) for nm = s(1)20 and r = 0(1)10. 


Table V. [np(1 — p)]! to 4D for m = 1(1)20 and p = .05(.05).5. 
Table VI. (pg)? to 4D for p = .005(.005).5. 
Table VII. e-™m?/x! for x = 0(1)---, to 4D for m = .1(.1)10(1)20. 


eo 
Table VIII. >> e-"m*/x! for the same ranges of the arguments as in 


z=mz! 


Table VII. 


t 
Table IX. y(t) = (2x)-te-*?, a/2 = f ¥(r)dr, and y™ (¢) to 4D or 4S 
0 


for ¢ = 0(.01)4(.05)5 and m = 1(1)6. 
p. 275. Selected important constants to 8D and their common logarithms 
to 7 or 8D. 


Table X(R). 1% and 5% points of the F distribution to 3 or more S for 
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the degrees of freedom 1(1)12, 14, 16, 20, 24, 30, 40, 50, 75, 100, 200, 500, 
for the larger estimated variance and 1(1)30(2)50(5)80, 100, 125, 150, 200, 
400, 1000, ~ for the smaller estimated variance. 

Table XI(R). 1%, 1% and 5% points of the z-distribution. (Z = 4 In F) 
to 4D for the degrees of freedom 1(1)6, 8, 12, 24, © for the larger estimated 
variance and 1(1)30, 60, «© for the smaller estimated variance. 

Table XII(R). The percentage points 100e of the Student-Fisher ¢ 
distribution to 3D for e = .001, .01, .02, .05, .1(.1).9 and the degrees of 
freedom 1(1)30, 40, 60, 120, ~. 

Table XIII. The cumulative Student-Fisher ¢ distribution to 3D for 
t = 0(.2)5 and the degrees of freedom 1(1)6, 8, 10, 15, 20. 

p. 285. A nomogram for 1 — (1 — p)*. 

Table XIV. (R)xe? satisfying P(x? = xo?) = eto 3D fore = .001, .01, .02, 
.05, .1, .2, .3, .5, .7, .8, .9, .95, .99 and the degrees of freedom 1(1)30. 

Table XV. e-* to 5D for x = 0(.01)3(.05)4(.1)6(.25)7(.5)10. 

Table XVI. m! to 5S and logio (m!) to 5D for m = 1(1)100. 

Table XVII. '(m) to 4D for nm = 1.01(.01)4.99. 

Table XVIII. logio '(m) to 4D for m = 1.01(.01)2. 

Table XIX. m! exact and (m!)—! to 5S for m = 1(1)20. 

Table XX. n? exact and n}, (102)? and 1000 n— to 5S for nm = 1(1)999. 

Table XXI. The natural trigonometric functions and the radian measure 
to 4D for the angle a = 0(1°)90°. 

Table XXII. In N to3D for N = 0(.01).99 and to 5D for N = 1(.01)10.09. 

Table XXIII. logio N to 4D for N = 100(1)999. 


c. © ¢. 


1212[K].—]J. H. CADWELL, “The distribution of quasi-ranges in samples 
from a normal population,” Annals Math. Stat., v. 24, 1953, p. 603-613. 


Let x1 S x2 S --- Sx, be the ordered values in a random sample of 
size m drawn from a normal population. The statistic w, = x,_, — Xp41 is 
called the range for r = 0, and a quasi-range for r = 1. The author gives an 
approximate, but apparently very accurate, method for evaluating the 
probability distribution function of w,. Of interest are the following results: 
(1) in section 3, the approximation for the p.d.f. of w, is given; (2) in section 
4, a comparison is made between the exact and approximate values of 
E(wo) and E(w) for n = 2, 8, 20, 30, 60, 100; (3) in section 5, the mean, 
variance, skewness, and flatness of the range are given for m = 20; (4) in 
section 7, there is a brief discussion of the efficiency of wo, wi, and we as 
estimators of the population standard deviation a; (5) in section 8, a table 
of values of moment constants to 3 or 4D and the .1, 1, 2.5, 5, 95, 97.5, 99, 
and 99.9 percentage points of w; are given for m = 10(1)30 to 2D; (6) in 
section 9, there is a brief consideration of the efficiency of estimates of o 
based on a linear combination of two quasi-ranges. The paper fills in a gap 
between published work for cases where m < 10 on the one hand and large 
on the other. 

BENJAMIN EPSTEIN 


Wayne University 
Detroit, Mich. 
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1213[K].—R. E. Ciark, ‘Percentage points of the incomplete beta func- 
tion,” Am. Stat. Assn., Jn., v. 48, 1953, p. 831-843. 


The table presented here gives to 4S the value of 


pb = P(N, X, a) 
defined by 
N (N 
ly = (*) oa — p)** = 1,(X, N — X + 1) 
for 


N = 10(1)50, X =1(1)N, a = .005, .01, .025, .05, 


where I,(X, Y) is KARL PEARSON’s incomplete beta function. Applications 
of the table are indicated. 

L. A. AROIAN 
Hughes Aircraft Co. 
Culver City, Calif. 


1214[K ].—W. J. Dixon, ‘Processing data for outliers,’ Biometrics, v. 9, 
1953, p. 74-89. 


Samples of size N are drawn from a mixture of normal distributions 
(1 — y)N(u, o?) + yN(u + Ao, o?) and (1 — y)N(u, 0?) + yN(u, d2o0*), re- 
spectively, in which y measures the contamination of the N(u, 0?) universe. 
Comparisons are made between the mean and median as estimates of yu; 
range and variance as estimates of dispersion. The bias and the MSE (mean 
squared error) for untreated (treatment refers to processing the data in a 
sample to remove outliers) samples (Table I), the ratios of the MSE of 
the mean and median of untreated data to the MSE of the mean for uncon- 
taminated data (Table I1), and minimum MSE of four procedures: use of 
X after treating for rejection at the significance levels, a = .01, .05, .10 and 
the use of median (Table III), serve as a basis for comparing the mean and 
median. The appropriate a to remove bias in s? and range (Tables IV, V, VI, 
VII) are used to compare range and variance. For most of the Tables 
N = 5,15; y = .01, .05, .10, .20, X = 0, 2, 3, 5, 7 with some modifications. 
The values known to be correct are those for \ = 0, the quantities for the 
mean in Table I, the results for no contamination in Tables II and III and 
the first line in Table V. Other results for N = 5, 15 are based on 100 and 
66 samples, respectively. The appendix gives critical values and criteria for 
testing for extreme values based on 17;; = (Xw — Xw-:)/(Xw — X41), 
a = .005, .01, .05, .10, .20, .30, N = 3(1)7 for rio, N = 8, 9, 10 for ri, 
N = 11, 12, 13 for ro1, N = 14(1)25 for ree to 3D. 

INGRAM OLKIN 

Michigan State College 
East Lansing, Mich. 


1215[K ].—AbBrauam Go vp, “Designing single-sampling inspection plans 
when the size is fixed,’’ Amer. Stat. Assn., Jn., v. 48, 1953, p. 278-288. 


Consider the sample size m fixed for a given single sampling plan. The 
problem is to choose an acceptance number s which yields the ‘‘best’’ pro- 
tection against misclassifying submitted lots. Two cases are considered: 
(1) single sampling-plans for placing a lot into one of two categories defined 
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by fraction defectives £; and 2, (2) single sampling-plans for placing a lot 
into one of m categories defined by fraction defectives 1, po, ---, Pm. In 
case (1), the “‘best’”’ plans are defined as those which minimize the sum of 
the producer’s and consumer’s risks. In case (2), the “‘best’’ plans are defined 
as those which minimize the sum of the probabilities of misclassifying 
the submitted lots. Given m categories defined by fraction defectives 
P1, D2, *-*, Pm, Tables 1 through 8 provide m — 1 acceptance numbers, 
c; (¢ = 1,2,---,m— 1) for single sampling-plans based on given fixed 
sample sizes 5(5)40, for acceptable lot qualities .01(.01).20, and for objec- 
tionable lot qualities .01(.01).40 respectively. The above tables should be 
very useful in quality control applications in which it is desired to divide 
the lot quality into more than two categories. 
G. W. McELrats 

Univ. of Minnesota 

Minneapolis, Minn. 


1216[K ].—E. L. Gras & I. R. Savace, “Tables of the expected value of 1/x 
for positive Bernoulli and Poisson variables,’’ Am. Stat. Assn., Jm., 
v. 49, 1954, p. 169-177. 


Let X be a random variable having the probability distribution 


@) Pix=8 =(%) peta - 2, 
where k = 1,2,---,n;q=1-—p;0<psl; 
(b) P{X =k} =e™(1 — e*)"m'/k!, 


where k = 1, 2, ---; m > 0. For case (a) Table I gives the values of 


n(n 
B(1/X|n,p) = (1 gE (j) eee 
to 5D for m = 2(1)20, p = .01, .05(.05).95, .99; m = 21(1)30, p = .01, 
.05(.05).50. The authors give an empirical approximation to (np — q)~ for 
certain values of the parameters which is accurate to at least 2.S. For case 
(b) Table II gives the values of 


E(1/X|m) = e™(1 — e*)> D m*k-/k! 
k=1 
to 5D for m = .01, .05, (.05)1(.1)2(.2)5(.5)7(1)10(2)20. 
INGRAM OLKIN 
Michigan State College 
East Lansing, Mich. 


1217[K ].—F. E. Grusss & H. J. Coon, “On setting test limits relative to 
specification limits,’ Industrial Quality Control, v. 10, No. 5, 1954, 
p. 15-20. 


This paper is concerned with the problem of setting test limits when 
specification limits are fixed and test measurements are subject to error. 
The authors consider test limits established in accordance with the following 
criteria: (1) producer’s and consumer’s risk are made equal; (2) the sum of 
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producer’s and consumer’s risk is minimized; and (3) the cost of incorrect 
decisions is minimized. As used here, consumer’s risk and producer’s risk do 
not have quite the same meaning that is usually associated with these terms 
in connection with sampling inspection. For the purposes of this paper, the 
consumer’s risk, designated as A, is defined as the joint probability of 
selecting a defective unit and judging it acceptable. Similarly the producer’s 
risk, designated as B, is defined as the joint probability of selecting an 
acceptable unit and judging it defective. With specification limits expressed 
as uw + ko, where a; is the product standard deviation, test limits are ex- 
pressed as uw + (ko, — bo.) where o, is the standard deviation of the errors 
of measurement. For criteria (1) and (2) above, tables of 6, A, B, and 
A + B are given to 4D with k = 1.5(.5)4, and r = .5(.25)1(1)10, where 
r = a,/o.. Examples are given which illustrate the use of these tables. 


A. C. CoHEN, Jr. 
University of Georgia 
Athens, Georgia 


1218[K].—J. M. Hammerstey, “On estimating restricted parameters,” 
Roy. Stat. Soc., Jn., s.B v. 12, 1950, p. 192-240. 


This paper contains a table of solutions of the following distribution 
problem: A philanthropist with u coins gives the first beggar he meets a 
number of coins which with equal probabilities may be 1, 2, ---, or uw. If 
the first beggar receives u — r coins, the next beggar will receive 1, 2, ---, 
or rf coins, each number being equally probable. If this is continued until 
the coins are exhausted, what is the probability that exactly » heggars will 
receive a gift? The probability P(u, v) satisfies the equation 


u—1 
P(u,vy) = YS w'P(r,v — 1) 
r=r—1 
with P(u, 1) = 1/u. P(u, v) is tabled for u and »v = 1(1)13. 
.. & &. 


1219(K ].—Hannes Hyrentus, ‘On the use of ranges, cross-ranges, and 
extremes in comparing small samples,”’ Am. Stat. Assn., Jn., v. 48, 1953, 
p. 534-545. 


A procedure is proposed using ranges, cross-ranges, and lower-extreme- 
differences as alternatives to the variance in testing variation homogeneity. 
The universe sampled is rectangular. Let sample 1 have N, items with lower 
extreme “,; and upper extreme v1, and let sample 2 have the corresponding 
No, U2, and v2. Designate the samples so that u; S wu. The ranges are 
Riu = 0; — & and Rey. = v2 — wu. The cross-range is Ro; = v2 — u;; and 
the lower-extreme-difference is Sx, = u2 — 4. 

Distributions of the test quotients T = S2:(Ry:), U = Re2(Ri1)— and 
V = Roi(Ri:)— are studied. With means and variances derived, tables are 
prepared showing upper 1, 5 and 10 percentage points of T, and showing 
upper .5, 2.5, and 5 and lower 95, 97.5 and 99.5 percentage points of U and 
V for Ni, N2 = 1(1)10 to 2D. 


T. A. BICKERSTAFF 


Univ. of Mississippi 
University, Miss. 
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1220(K ].—R. Latscna, “Tests of significance in a 2 X 2 contingency table: 
Extension of Finney’s table,’”’ Biometrika, v. 40, 1953, p. 74-86. 


Consider a 2 X 2 contingency table arranged so that the sum of the two 
values in the first row is A and the sum of the two values in the second row 
is B where A = B. Let a be value in the first row and first column while } 
is the value in the second row and first column. For given A, B, a, 6, the 
table presented is entered in the section for A, the sub-section for B and 
the line for a; then the main body of the table shows in bold type the 
appropriate significance points for b. If the observed value b is equal to or 
less than the bold-faced integer in the column headed .05, .025, .01 or .005 
(only values considered), then a/A is significantly greater than b/B (single- 
tail test) at these probability levels. This same condition on 6 furnishes a 
two-tail test of whether a/A differs from b/B at twice the probability value 
heading the column. A dash or absence of entry for specified A, B, a indi- 
cates that no 2 X 2 table in that class can show a significant effect at that 
level. The true significance level is at most equal to the probability value 
heading the column and in some cases may be noticeably smaller. For each 
table entry considered, the true probability value is given to 3D in small 
type next to the bold-faced integer with which 5 is compared. The table 
presented in this article covers the cases A = 16(1)20, B = 2(1)A, A 2a 
= (smallest integer not yielding all dashes). Previously FrnNEy! obtained 
a similar table for A = 9(1)15. 

J. E. WALsH 
U. S. Naval Ordnance Test Station, Inyokern 
China Lake, California 


1D. J. Fryney, “The Fisher-Yates test of significance in 2 X 2 contingency tables,” 
Biometrika, v. 35, 1948, p. 145-156 [MTAC, v. 3, p. 359]. 


1221[K ].—M. Masuyama & Y. Kurorwa, “Table for the likelihood solutions 
of gamma distribution and its medical applications,’’ Union of Japanese 
Scientists and Engineers, Reports of Statistical Application Research, 
v. 1, 1951, p. 18-23. 


For the I'-type frequency function, 
fy) = (y/a)?— exp (— y/a) dy[al'(p)F" if y20 
= \0 if y <0 


the maximum likelihood estimates, 6 and 4, of the two parameters, p and a, 
are the solutions of the equations, 


In p — ¥(p) = In (A/G) = g() 
4=A/p 


in which y is the digamma function and A and G are the arithmetic and geo- 
metric means of the sample respectively. To facilitate the approximate 
solution of the equations, the authors have tabulated x log (A/G), g(x), 
log (A/G) and A/G to 8S or 7D for x = .1(.05)3(.1)5(.5)10, 15, 20(10)50. 


No@ , 
For large samples of N, — and No;? are given to 8D for the same values of x. 


G.ac 
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1222[K].—NBSCL Probability Tables for the Analysis of Extreme-Value 
Data. NBS Applied Math. Ser. No. 22. U. S. Gov. Printing Office, 
Washington, D. C., 1953. iii + 32 p., 22 XK 20 cm. Price $0.25. 


There are many engineering and scientific problems in which the random 
variable being observed is a largest (smallest) or more generally m’th largest 
(m’th smallest) value. The 15 page introduction by E. J. GuMBEL, who 
pioneered in the field of extreme values, is a succinct and excellent survey of 
current knowledge in the field. In the first 9 pages of this introduction he 
sketches the theory and carries out the analysis of a problem in detail. The 
remainder of the introduction is devoted to listing the functions tabulated, 
giving a history of the genesis of the tables, and describing methods of 
computation. A useful bibliography completes the introduction. 

If one carries out a suitable transformation on the largest value in a 
sample one gets a random variable which is called the reduced largest value. 
The asymptotic distribution of this reduced largest value is (for a wide class 
of common distributions such as the exponential, normal, and chi-square 
distributions) given by the cumulative distribution function 6, = exp (—e~”) 
and associated probability density function g, = %,’ = exp (— y — e~”). 
Table 1 gives 4, and g, for y = — 3(.1) — 2.4(.05)0(.1)4(.2)8(.5)17 to 7D. 
Table 2 gives the inverse of the cumulative probability function of ex- 
tremes y = — log, (— log, ,) for &, = .0001(.0001).005(.001).988(.0001) 
-9994 (.00001).99999 to 5D. Table 3 gives the probability density function ¢, 
as a function of the cumulative distribution function ®,[ ¢, = — , log #, ] 
for , = 0(.0001).01(.001).999 to 5D. The reduced m’th largest value 

Ym ™ 

has associated cumulative distribution function 4, = f ai Di 
X exp (— my — me~)dy. Table 4 gives probability points y,, for the reduced 
m’th largest value, Ym = Ym(Pm), for m = 1(1)15(5)50; & = .005, .01, .025, 
.05, .1, .25, .5, .75, .9, .95, .975, .99, .995 to 5D. Tables 5 and 6 deal with an 
appropriately reduced range R. The cumulative probability distribution 
function ¥(R) is given by ¥(R) = 2e-#*K,(2e-?*) and the probability den- 
sity function ¥(R) = W’(R) is given by ¥(R) = 2e-®Ko(2e**). Ko and K, 
are the modified Bessel functions of the second kind of orders 0 and 1, respec- 
tively. Table 5 gives ¥(R) and ¥(R) for R = — 4.6(.1) — 3.3(.05)11(.5)20 
to 7D. In Table 6, the reduced range R is considered as a function of Y, i.e., 
as the solution of the equation VY = 2e-#*K,(2e!*). This is given in Table 6 
for Y = .0001(.0001).001 (.001).01 (.01).95 (.001).999 (.0001).9999 to 4D, ex- 
cept for ¥ = .0001 and = .999 to 3D. 

The tables are the fruits of the painstaking labors of many people. 
Chief among these are E. J. GUMBEL who took part in the detailed planning, 
J. ARTHUR GREENWOOD, JULIUS LIEBLEIN, and H. E. Sauzer. The prepara- 
tion of the tables in their final form was carried out under SALZER’s direction 
at the Computation Laboratory of the National Bureau of Standards. 
LIEBLEIN assisted in the preparation of the Introduction. 

BENJAMIN EPSTEIN 


Wayne University 
Detroit, Mich. 
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1223[K }.—NBSCL Tables to Facilitate Sequential t-Tests. NBS Appl. Math. 
Ser. No. 7. U. S. Gov. Printing Office, Washington, 1951. xix + 52 p. 
Price $0.45. 


The sequential ¢-test is designed to test the hypothesis that the mean of 
a normally distributed variate with unknown mean and unknown variance 
has a given value. The present tables are designed for use with a probability 
ratio test based on the likelihood ratio for the non-central ¢ to the central ¢, 
which thus differs from the one originally proposed by WALD.' However by 
a simple adjustment of the parameters these tables can also be used for 
Wald’s test. 

The test criterion after the n-th observation is Z, = [x (xa — lo) P 


4 ‘>> (%a — Io)? }-! in which Jp is the value of the mean iontiin by the 


a=l 
null hypothesis. The critical values of Z, are solutions of the equation 
&Z & 
L=In P( 1/2, > ) = + in which F(a, 6b; x) is the confluent hyper- 
geometric function, 


2 (6) (a +f) x 
a (aE +f) ji!’ 


§ is the amount in standard units by which the mean under the alternate 
hypothesis differs from J) and L is either In {(1 — 8)/a} or In {8/(1 — a)} 
in which a and £ are the risks of errors of the first and second kinds respec- 
tively. Values of Z, are given to 1D for 6 = .1, to 2D for 6 = .2(.1).5 (to 
3D for L < 0 and 6 = .2) and to 3D for 6 = .6(.1)1(.2)2, 2.5 for + L = 2, 
In 19, 3, 4, In 99, 5, 6, 7 and m = 1(1) S 200. (Values of are carried as far 
in each case as is judged likely to be useful in practice but never beyond 
200.) There are some smaller auxiliary tables. One of these gives the values 
of a and 8 to 3S for the L’s listed ; others are to assist in reaching a decision 
if ~ goes beyond the values in the table; and still others give approximate 
upper bounds for the error in linear interpolation with respect to L or S. 


cc &."<. 
1A. WALD, Sequential Analysis. New York, 1947, p. 204-207. 


1224[K ].—B. M. SEELBINDER, “On Stein’s two-stage sampling scheme,” 
Annals Math. Stat., v. 24, 1953, p. 640-649. 


Stern! has proposed a two-stage plan for sampling from a normal popula- 
tion of unknown variance o? in order to estimate the expectation with a 
confidence interval of preassigned length 2d and confidence coefficient 1 — a. 
The experimenter sets the size , of the first sample, while the size n — m,; 
of the second sample depends on the variability observed in the first. The 
present paper is concerned with the problem of choosing m; with an eye to 
minimizing E(n). SEELBINDER tables E(m) for some of the following values: 
a = .1, .05, .02, .01; d/o = .01(.01).1(.1)1; m1 — 1 = 5, 10(10)60, 80, 120, 
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240. An approximation based on the normal distribution is found to be 
adequate. 
J. L. Hopcgs, Jr. 


Univ. of California 
Berkeley, Calif. 


1 CHARLES STEIN, “A two-sample test for a linear hypothesis whose power is independent 
of the variance,” Annals Math. Stat., v. 16, 1945, p. 243-258. 


1225[K,N ].—H. Sremnnaus, Table of shuffled four-digit numbers. Rozprawy 
Matematyczne, No. 6, Warsaw 1954, 46 p. The introduction is in Polish, 
Russian and English. 


This table of 10,000 random 4-digit numbers differs from the usual 
random number lists in that no two 4-digit numbers are equal. In other 
words the table, when read in the usual order, gives one of the 10,000! 
permutations of the numbers 0000-9999. This feature allows the list to be 
used in problems in which samples are drawn without replacement. One 
may also read the digits vertically by fours and obtain samples which are 
drawn with replacement. 

The table was produced by hand. To begin with a table of 10,000 four- 
digit numbers U,(m = 0(1)9999) was produced by 


U, = 4567U,_1 (mod 10,000) 


with U» = 0000. This table was arranged in a square 100 X 100 and then 
subjected to a number of randomizing transformations fully described in 
the introduction. Whether the 25 pages which finally resulted pass any 
of the standard tests of randomness is not indicated. Apparently the author 
has made no tests whatever. The reviewer subjected the first page to a 
frequency test which gave a x? with a probability of only 2 percent. 
The list could be useful for ordering the retirement of bonds. 
D. H. L. 


1226[K ].—G. M. Tuompson, “Scale factors and degrees of freedom for 
small sample sizes for x approximation to the range,’”’ Biometrika, v. 40, 
1953, p. 449-450. 


Consider m independent samples of size m from a normal population, 
mean yw, variance o*. Let w be the range of a sample, @m,, the mean range. 
PaTNAIK! found the approximation 

@m,n = cxv— to 
where c¢ is a scale factor, vy the degrees of freedom for x. Evidently Gm, n/c 
is equivalent to the usual standard deviation estimator with v degrees of 
freedom. Thompson tables c and »v to 4D for m = 1, m = 2(1)10. He also 
gives the 95% confidence limits on uw using w, ie. X + k.osw to 4D, for 
n = 2(1)10. Agreement with results obtained by the use of Student’s ¢ is 
excellent. 

H. A. FREEMAN 

Mass. Inst. of Technology 
Cambridge, Mass. 


1P. B. Patnalk, “The use of mean range as an estimator of variance in statistical tests,” 
Biometrika, v. 37, 1950, p. 78-87. 
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1227[L ].—F. DrMaceio, A. Gomza, W. E. Toomas & M. G. SALvaport, 
“Instabilita laterale di travi inflesse e compresse,’’ Accad. Naz. dei Lincei, 
Atti Rend., s. 8, v. 12, 1952, p. 524-529. 


This paper summarizes the results of a number of technical reports 
published by the Dept. of Civil Engineering, Columbia University. The 
critical loads for the lateral buckling of several types of beams and loadings 
are given. These are the lowest eigenvalues of 5V = 0, where 


l l I 
V =3B f (u’’)*dz 4- pie f (B’’ "dz + an f (p’)*dz 
0 0 0 
I I 
+ f MBy"dz — 3P f (u’)*dz. 
0 0 


B, D, h, H and P are constants. Table 1 gives the results for I-beams under 
unequal end moments, tabulates eigenvalues K = M.//(BC)! to 2D asa 
function of P/a? = 2PC/hD = .1 (irreg.) © and r = M,/M; = — 1(.1)1, 
where H = C — PI,/A, M(s) = M.[r+ (1 — r)z/l], Mi = rMse, u(0) 
= p’"(0) = n()) = vw" (1) = 0, B(0) = 6’ (0) = BY) = B”(D) = 0. Table 2 
gives the solution of rectangular beams under bending and compression, 
corresponding to D = 0, u(0) = u(J) = w’’ (0) = »” (D) = 0, B(0) = BY) = 0. 
K = M.I/(BH)}' given to 2D for p = PPx?/B = 0(.2)1; and r = M,/M;2 
= — 1(.1)1. Table 3 gives solution to rectangular cantilever beam with 
shear and axial load at the free end corresponding to D = 0; r = 0; M(z) 
= Qz = Mz2/land w(l) = w’(l) = vw" (0) = Bu’ (0) + Pu’(0) + QB(0) = 0, 
B(l) = B’(0) = 0, g = QP/(4.013(BH)') given to 4D for p = 4PPr?/B 
= 1(— .2) — .4, — 1. Solutions of this problem (originated by PRANDTL 
and MICHELL) for many other end conditions have been discussed in 
TIMOSHENKO.! 
Y. C. FunG 


Calif. Inst. of Technology 
Pasadena, Calif. 


1S. TimoOsHENKO, Theory of Elastic Stability. New York, 1936. 


1228[L ].—Stic Exe .6r, ‘“‘Theory of electromagnetically delayed telephone 
relays,” Ericson Technics, v. 9, 1953, p. 141-224. 


On p. 221-224 there is a 4D table of the functions 


sinhx 3coshx 3sinhx 2(coshx — 1) 











cosh x, 7% - lee ge x , 
12sinhx  24(coshx — 1) 
= 4 x 


for x = 0(.01)1-5. Graphs of these functions (in the same range) are given 
on p. 161. 


A. E. 
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1229[ L ].—Ernst GLowatTZkI, ‘‘Tafel der Jacobischen elliptischen Funktion 
@ = am (mK/n),”’ Bayer. Akad. Wiss., math.-naturw. Kl., Abhand- 
lungen, Neue Folge, heft 61, 1953, 27 p. 
3D tables of am(u, k) in degrees for u = mK/n, k = sin 0, nm = 2(1)12, 

nm = 1(1)m — 1, 6° = 0(1)90. First differences are also given. The entries 

were computed from Legendre’s tables, and were checked against the 

Smithsonian Elliptic Function Tables (see MTAC, v. 3, 1948-49, p. 89, 

RMT 485). 

A. E. 


1230[L].—Louts Rosin & ALFREDO PEREIRA-GoMES, “L’antenne bi- 
conique, symmétrique, d’angle quelconque,”’ Ann. d. Télécommunications, 
v. 8, 1953, p. 382-390. 


The authors discuss the equations 





P,(cos 0) + P,(— cos@) = 0 or F(- ‘5 £14: costo) = 0 
n 


1 — n 


#42 scot) = 0 





P,(cos 6) — P,(— cos@) = 0 or (45 


where P,, is the Legendre function, F the hypergeometric series, @ is given, 
and n is to be determined. 

On p. 390 they give tables (mostly to 3D) of the first nine roots of the 
first equation, and of the first eight roots of the second equation, when 
6 = 2/12, x/6, 2/4, 2/3, 5/12. 

A. E. 


1231[L ].—M. Rotuman, “The problem of an infinite plate under an inclined 
loading, with tables of the integrals of Ai(+«) and Bi(+<),’’ Quart. Jn. 
Mech. Appl. Math., v. 7, 1954, p. 1-7. 


For the definition of the Airy integrals Ai(x) and Bi(x) see MTAC, v. 1, 
1944, p. 236. 


Table 1. 7D values of f Ai(é)dt for x = 0(.1)7.5 with 6,2. 7D values 
0 
of f Bi(é)dt for x = 0(.1)2 with 6,.2 and 7‘. 
0 
Table 2. 7D values of f Ai(—dé)dt for x = 0(.1)10 with 6,,2 and 4. 
0 


Table 3. 7D values of f Bi(—é)dt for x = 0(.1)10 with 6,,2 and y‘. 
0 
A. E. 


1232[L].—F. G. Tricomr, “On the statistical distribution of mutant 
bacteria,’’ Bull. Mathem. Biophysics, v. 15, 1953, p. 277-292. 


The author gives 5D tables of 


“ " dG 
G)=E eae CO=G. GO = [cman 
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for = — 2(.1)2. Graphs are given of G(#), and also of the curve represented 
parametrically by 


‘“ i 
x=G(), y= [ ea | exp {alG*(t) — 1G(4)]]} 
for a = 1 anda = .01. 
A. E. 


1233[L ].—E. F. M. vAN DER HELD, “The contribution of radiation to the 
conduction of heat. II. Boundary conditions,” Appl. Sci. Research A, 
v. 4, 1953, p. 77-99. 


The Appendix (p. 92-99) contains a collection of formulas and some 
numerical tables concerning the function 


K,(x) = f t-"e-* "dt 


which is closely related to the complementary incomplete gamma function. 
Table VII gives values of Ki, Ke, 2K3, 3K4, 4K, to 4D or 4S for 
x = 0(.01).2(.02)2(.5)5(1)8. 
Table VIII gives values to 2 to 4S of 


f t—1e-** log t dt, f K,(x — )Ke(é)dt, nm = 1(1)5 
1 0 


for x = 0, .01, .05, .1, .25(.25)1(.5)5(1)8. 
Table IX gives values to 2 to 4S of 


f K,(D — x + )K2(0)dt, ‘ K(x +)Ks()dt, D=2,4, n= 2(1)5 


in varying ranges, and at varying intervals, of x. 

V. Kourcanorr (see MTAC, v. 3, p. 307, RMT 569) has tabulated 
Ke, 2K3, $ — 2K, $x + 4K; for x = 0(.02)2. Earlier tables of K,(x) are 
listed in FMR Index, sec. 14.83. 

A. E. 


1234[S ].—Otro EMERSLEBEN, “Das elektrostatische Selbstpotential aqui- 
distanter Ladungen auf einer Kreislinie,’’ Math. Nachr., v. 10, 1953, 
p. 135-167. 


If 2” charged particles are equally spaced on a circle of radius r with a 
charge (—1)* at Pi, 1 < k < 2m, then the self-energy of the system is given 
by the finite sum 


n (n) 
sb =" (—)hesc (buk/n), 
k=1 
where (m) indicates that the last term must be multiplied by }. In a unit lat- 
tice » = mr and the power series for the cosecant leads to the decomposition 
2 - aPo + aP1, 


where ,®o comes from the first term of the power series and is similar to 
the potential of a linear configuration. Both expressions ,%9 and ,®, involve 
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(n) 
the sums }> (—)*k-* for various integer values of s and their asymptotic 
r=1 


behavior for large m is studied. The term ,®; arises from the curvature of 
the circle and is O(n-") as n > ~, while ,B) = — (2 log 2)n + O(n). The 
O-term of ,®o is also analyzed in somewhat more detail and a table of > 
is given for m = 1(1)12. 

Similar results are obtained for the problem of equally charged and 
equally spaced particles on a circle. The self-energy in this case is 


e ae = 
a? = b au csc (xk/n), 


and a table is given for » = 1(1)30 when nm = 2zr. 

The alternating sums > {—1)*k™ reduce to polynomials in with ra- 
tional coefficients when a” % an integer. These are given explicitly for 
m = 0(1)7. Finally, the graphs of the seven sums > (—1)*k*, m = 1(1)7, 


are plotted as functions of s for — 2.6 < s < 3. 


Tom M. ApostTo. 
Calif. Inst. of Technology 
Pasadena, Calif. 
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In this issue references have been made to errata in RMT 1205. 


239.—NBSCL, Tables of 10*. (Antilogarithms to the Base 10.) Applied Math. 
Ser. No. 27. Washington, 1953. 


In working with 10-figure logarithms recently, I noticed a discrepancy 
between a value given in this new table of 10", and a value from GUILLEMIN.! 
Since Guillemin gives 13 figures for 107 for x = 0(.0001).6999, and 12 figures 
for x = .6999(.0001).9999, and DEpREz? gives 14 places 1(1)9999, every 
tenth value in v. 27 may easily be checked directly from these tables. 

In checking every tenth value for the block from x = .40000 to .50000, 
I find 128 values in error in the last figure; 127 values should have the last 
figure raised by one, and one value, that for .49270, should have the last 
figure reduced by one. 

It seems rather surprising that a corps of experienced computers, working 
with the most modern calculating instruments, should spend the amount of 
time and labor represented by this table simply to smooth out a 200-year-old 
table without even checking the accuracy of that table. Only the expenditure 
of a very little extra time would have been necessary to use the Sundstrand 
machine to subtabulate Deprez’s table to tenths and thus to get a 14-place 
table of 10°, which, if rounded to ten places would yield a table far superior 
to Dodson. Reference is made in the introduction of this volume to the 
Deprez table, so it is certainly known and available to the personnel of the 
laboratory. 

In line with this reasoning, by using the Deprez values for the block 
-4680 to .4700 in conjunction with Comrir’s simple, accurate, and very 
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powerful system of end-figure interpolation to tenths (Supplement to Nautical 
Almanac, London, 1931), together with an old Ellis two-register machine, 
I reproduced the block .46800 to .47000 (p. 236) to 14 places. Rounding this 
to ten places and comparing page 236 with these values, I find 23 further 
errors in the 200 entries, besides the four errors already found in this block. 
Judging from this sample and from the test above, it would appear that 
about 13 per cent of the entries in Table I are in error by being one too low 
in the last digit. 

This same volume contains a very fine table, Table II, which gives 16- 
place values for 107 for x = 0(.001).999. Every 100th value in Table I may 
be checked easily by direct comparison with the second column of Table II. 
In the block .40000 to .50000 there are eleven of this type of entry in Table I 
in error, as may be directly ascertained by checking in Table II. 

It is to be hoped that in the near future the Laboratory will issue a cor- 
rected version of this Table. 

A list of the errors referred to above is appended herewith. 

The final digit in 107 is too low by a unit for the following values of x. 


40020 41690 43920 46070 46916 48190 
40100 780 980 080 921 260 
40190 830 990 180 928 270 
320 900 44010 240 930 380 
390 930 090 700 936 390 
440 950 140 710 938 520 
450 42110 180 750 944 570 
570 160 300 780 963 670 
590 240 400 800 965 730 
620 260 450 804 990 750 
700 520 540 825 47100 810 
790 740 610 832 160 910 
840 810 620 834 170 920 
41220 900 680 837 190 49050 
300 43020 690 838 250 180 
350 050 740 852 550 530 
500 070 780 855 580 650 
510 160 800 860 640 750 
520 350 45110 863 750 780 
530 480 190 867 790 840 
550 530 260 868 850 870 
570 670 270 903 870 880 
590 730 360 908 940 930 
620 840 490 909 950 980 
640 910 46030 912 48140 990 


FREDERICK W. HOFFMAN 
689 East Avenue 
Pawtucket, R. I. 


1A. GuILLEMIN, Tables de Logarithmes @ 3 quatrades et Nombres Correspondants avec 
12-13 chiffres. Paris, 1912. 

2F, Deprez, Tables for Calculating, by Machine, Logarithms to 13 Places of Decimals. 
Berne, 1939. 
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188(D,F].—D. H. Lenmer, Roots of Unity and Cyclotomic Periods. Tabu- 
lated from punched cards and deposited in the UMT Fixe. 


In UMT 145 [MTAC, v. 6, p. 102] we gave values of the cyclotomic 
cosines 2 cos 2rk/p for k = 1(1)p — 1, for all primes p < 97 to 20D. We 
now have a companion table of the cyclotomic sines for the same values of 
k and p to 12D, and therefore a table of the exponential exp (2rik/p) for 
these values of k and # or in other words of the primitive p-th roots of unity 
for p a prime <100 to 12D. 

If we write p = ef + 1, then the roots of unity may be summed in e 
“periods” of f roots each according as the index of k with respect to a 
primitive root g [these primitive roots were chosen as in JACOBI" ] is con- 
gruent to 0, 1, 2, ---,e — 1 (mod e). These periods n, (s = 0,1, ---,e — 1) 
are roots of period equations of degree e with integer coefficients, whose 
coefficients, discriminant and sum of whose powers are of interest in the 
theory of numbers. The present table tabulates 7, for s = O(1)e — 1, 
e = 5(1)16, 18(2)26, 32 for all primes of the form p = ef + 1 < 100. 


im. 
1C. G. J. Jacost, Canon Arithmeticus. Berlin, 1839. 


189[F].—A. Ferrier, Etude de 5n? — 1 et 5n* — 1. 35 typewritten pages on 
deposit in UMT Fite. 
This study is in two parts. The first part gives the solutions m of the 
congruences 
Sn? 


1 (mod P) P= 10x +1 < 11,370 
and 
5n4 


1 (mod P) P = 10x +1 < 10,000 


The second gives the factorizations of 5? — 1 and 5m‘ — 1 for n < 5000 
and m < 1000 respectively. All primes g above 107 which appear have been 
checked by showing that g divides 2" — 2. This was done on the EDSAC 
with the cooperation of J. C. P. MILLER. 

A. FERRIER 
Collége de Cusset 
Allier, France 


190[F].—A. GLopEn, Factorization of N* + 1 for isolated values of N between 
500 and 4000. Three typewritten pages on deposit in the UMT FILe. 


This table is an extension of a table of CUNNINGHAM! for N < 200 and 
a table by the author for N < 500 [MTAC, v. 7, p. 34, UMT 153]. In the 
present table 53 complete factorizations of N* + 1 which are all new with 
the exception of N = 512 and N = 1024. 
A. GLODEN 
11 Rue Jean Jaurés 
Luxembourg 


1A. J. C. CunninGHAM, Binomial Factorisations, v. 1. London 1923, p. 140-141. 
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Edited by the Staff of the National Applied Mathematics Laboratory of the National 
Bureau of Standards. Correspondence regarding the Section should be directed to Dr. E. W. 
CANNON, 415 South Building, National Bureau of Standards, Washington 25, D. C. 


TECHNICAL DEVELOPMENTS 


A Punched Card Machine Especially Equipped 
for Scientific Computations 


The author believes that a description of the punched-card machine now 
in use at the Norwegian Defence Research Establishment might be of some 
general interest. 

The nucleus of the installation is a type C3 multiplying punched-card 
machine, purchased from the French Compagnie des Machines Bull. The 
original plugboard supplied at the time of purchase has been replaced by 
another almost double the size, and, by making certain other additions, the 
range of usefulness of the machine has been extended to include many types 
of calculations which are not normally practicable with punched-card 
machines. 

The machine has 96 counting wheels, which are connected in groups of 
three by carry-over. Connection of the groups is determined by the plugging 
of the carry-over. Four 3-digit counters may be used for a multiplier and 
four others for a multiplicand, thus allowing two 12-digit numbers to be 
multiplied. The product is formed by the usual left and right hand principle, 
the two partial products being finally added together. It will be seen that 
even when two 12-digit numbers are multiplied, there are eight 3-digit 
counters which are not used. Although the counters do not permit subtrac- 
tions to be performed directly, a number may be subtracted by adding its 
nines complement, and using ‘‘end-around” carry. 

The operations of multiplication, addition, subtraction, conversion of a 
number into its nines complement and clearing may be performed. In any 
cycle only one type of operation can be performed. In order to transfer a 
number from one counter group to another, 48 change-over relays are pro- 
vided, each relay having three sets of contacts. 

A programme is limited by the 44 cycles which are built into the ma- 
chine. However, the cycles may be connected in any desired sequence and 
used over and over again. With each cycle is associated a total of 4 “‘com- 
mand” hubs, which are connected to the positive side of the 48-volt supply 
during the time the cycle is in progress. The operation to be performed dur- 
ing a cycle is determined by plugging from the hubs, the change-over relays 
being energized in the same way. 

It is of interest that a multiplication must be completely programmed, 
one cycle for each digit in the multiplier and one for the addition of the 
partial products. Two command hubs are used for each digit, one determin- 
ing the digit selected from the multiplier and the other determining the 
amount of left shift of the partial product. Furthermore, any digit or digits 
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of the multiplicand may be suppressed during multiplication with any digit 
of the multiplier. This feature permits several trick-pluggings, e.g., several 
products may be formed at the same time, and it may even be possible to 
obtain the sum of several products in one operation. 

In the event that more than one multiplication occurs during a problem, 
a uniselector is provided. The uniselector is stepped during each multiplica- 
tion, usually in the cycle where the partial products are added. The impulse 
which starts the cycle following a multiplication is led to the input of the 
uniselector, and the cycle which is required to follow the multiplication is 
then initiated from the output hubs of the uniselector. These hubs, 50 in 
number, may be used for a variety of other purposes as well. 

Eight sign-control units have been connected to eight experimentally 
selected counting wheels, those which experience in use of the machine 
have appeared to be most useful for this purpose. Each unit has a change- 
over relay, six relays having two sets of contacts, and the remaining two 
having three sets of contacts. The sign control takes up one position if the 
associated counting wheel indicates a nine and assumes the alternative 
position when any other number is indicated. 

By means of a comparison unit, provision is made for comparison at the 
upper and lower brushes of two 6-digit numbers from two consecutive cards. 
Mention must also be made of ten relay units which can remember a pulse. 
Six of these units are cleared at the beginning of a card cycle, and the other 
four by impulses sent to the clearing hubs either by x punches or in any 
cycle by pulses from the command hubs. In addition, there have been in- 
stalled eight units similar to the first six, except that they clear in the card 
cycle when the zeros are passing the brushes. If these units are set up from 
the lower brushes they work exactly as the other six, but by setting them up 
from the upper brushes other obvious effects are obtained. The outputs from 
the units consist of a positive voltage on an output hub which is normally 
connected to the bobbin of one of the change-over relays. The sign-control 
relays, and the relays controlled by memory units, may be used to switch 
the succession of cycles. Thus, although the machine is not directly designed 
for division, this operation can easily be programmed. 

Numbers may be extracted from the machine in three different ways. 
They can be punched into the cards supplied to the machine, or into new 
cards passing through an auxiliary punch, or they can be recorded by 
teleprinter. In the first instance the machine must stop, but this may be 
avoided in the other two cases by inserting a 24-digit relay-storage unit. 

In order to check the complicated plugging required in some problems, 
a panel with 300 indicator lamps, to show the state of all important relays, 
has been added. A special feature permits the running of single cycles, 
allowing observation of the behaviour of the machine during this period. 

Further technical details are likely to be of interest only to those actually 
using the equipment, and the author presents now an outline of some 
problems solved. 

(a) Calculation of primes. This problem was attempted only as a demon- 
stration test. The method of successive division was used. Observation was 
made after each division to see if the remainder was zero or not. Division 
in each case was performed by subtracting the divisor, 10 times the divisor, 
100 times the divisor, etc., until the dividend became negative; division 
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proper was then started from this point. In this way a minimum number of 
cycles was used. The machine ran continuously for 30 hours, and produced 
the primes up to 3600 on the teleprinter output. A single punched card was 
used, giving 7 as the first number to be tested. 

(b) Calculation of ballistic curves. The ballistic curves for a slow rocket 
projectile were obtained from 


= — ket 

7 = — km —g 
g=¢ 

y= 

v= 2+ 7’. 


In the case considered kv = a + b(# + 7”) and the simple integration 
formula 

£1 = €1 + 2héo 
could be used. 

The initial values of the first two instants and the time of flight were 
punched into two cards for each curve. The values of t, , 7, x and y at second 
intervals during the total time of flight of 20 seconds were given in tabular 
form by the teleprinter. Eighty curves were computed in a little over 
six hours. 

(c) Solution of linear equations. The solution was obtained by Verzuh’s 
method, (see MTAC, v. 3, p. 453-462), but, because of the extra facilities 
available, the process was considerably simplified. 

The first equation was divided by its leading coefficient, and two cards 
were punched by the auxiliary punch for each coefficient. By means of a 
sorter these cards were then mixed with the cards of the other equations, 
and one card was punched for each of the coefficients in the resulting equa- 
tions, where one unknown was eliminated. The reason for computing two 
cards in the first case was that, whilst one set of cards was used by the 
machine for eliminating an unknown from one equation, the other set could 
be mixed with the cards giving the next unknown. Normally several equa- 
tions are solved simultaneously. Whenever the machine has finished the 
calculations for one equation, the teleprinter prints the sum of the cards 
just punched. As there is an auxiliary card for each equation, containing 
minus the sum of all the other coefficients in the equation, the teleprinter 
would always print 0 or —1 if no round-off errors existed. Each equation 
was preceded by an X-punched master card, different for the first equation 
and the other equations. This card not only caused the teleprinter to print 
the contents of the counter giving the sum of the previously punched cards, 
but also cleared this counter and switched the machine to the correct pro- 
gramme. No manual switching was therefore required. Ten sets of four 
equations with four unknowns were solved in two hours by this method. 

The author believes that it would not be worth while to construct a more 
flexible electromechanical computor than the one described, because of the 
inherent slow speed of such punched-card machines. The introduction of 
electronic techniques into the conventional punched-card system is a stage 
which will probably be superseded very soon by the use of a true all-purpose 
electronic computor. Such a machine may, if so desired, have punched cards 
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as input and output media. For this reason the author has not incorporated 
an electronic multiplier in the machine described. The cost of the equipment 
is about 20,000 dollars, and it is already possible to buy an electronic general- 
purpose computor for about this price. It would seem therefore that the 
era in which punched-card machines have been applied to scientific com- 
putations will very soon be past. 
Jan V. GARWICK 

Norwegian Defence 

Research Establishment 
Kjeller, Norway 


BIBLIOGRAPHY OF CODING PROCEDURE 


The material described below is among that which has been added to the 
collection at the National Bureau of Standards Computation Laboratory. 
A similar collection is available at the Numerical Analysis Research, Univ. 
of California, Los Angeles 24. 

Material for inclusion in these coliections should be sent to the Computa- 
tion Laboratory, National Bureau of Standards, Washington 25, D. C., for 
the attention of J. H. WEGsTEIN. 


56. J. W. Backus, ‘““The IBM 701 Speedcoding system,’’ Assoc. Computing 
Machinery, Jn., v. 1, 1954, p. 4-6. (Cf. item 58 of this Bibliography.) 

57. H. H. GoLpDsTINE, ‘Some experience in coding and programming with 
the Institute Computer,”” Argonne National Laboratory, Proceedings of 
a symposium on large scale digital computing machines, August 3-5, 1953, 
ANL Report 5181, p. 273-278. 

It is pointed out that it is not possible, with present day equipment, 
to divorce the engineering and mathematical difficulties in program 
preparation. Limitations on read-around ratio imply restrictions on 
coding ; on the other hand, mathematically redundant checks are useful 
both in code checking and problem running. 

A code-checking scheme recommended is first the checking of the 
arithmetics of the code (e.g. sub-routines purporting to represent a 
function accurately or approximately) and then the logics of the code. 
In the first step a specialization of parameters is often helpful (too great 
specialization must be avoided—cf. R. A. Brooker, S. GrLi, D. J. 
WHEELER, MTAC, v. 6, p. 112-113). In the second stage it is useful to 
record the path of the control through the sub-routine complex. Fast 
printout of the full memory, to be contemplated at leisure, away from 
the console, has proved advantageous. 

The importance of keeping an adequate diary describing the day-to- 
day development of the solution of a problem is stressed. 

It is interesting to note that it has been decided to add hardware to 
the machine to detect overflows or spills. 

By 
NBSCL 
58. INTERNATIONAL BusINESS MACHINES Corp., Bibliography on the Use of 
IBM Machines in Science, Statistics, and Education. New York, Janu- 
ary 1954. 
59. W. L. VAN DER POEL, “‘Dead programmes for a magnetic drum automatic 
computer,” Applied Sci. Research, v. 3B, 1953, p. 190-198. 
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Dead programming here refers to programs stored in permanent 
memory and this paper points out the value of such programming. 


. C. Ross, C. S. FLuKE, H. S. Toney, Handbook for OARAC Programming 


and Coding. Wright Air Development Center, Nov. 1952. 

This is a concise set of notes for guidance and reference in pro- 
gramming for the OARAC. It includes the order code, a glossary of 
terms, and examples of simple codes. 

Programming Univac fac-tronic Systems. Instruction Manual 1. Rem- 
ington Rand Inc., Philadelphia, 1953, 249 pages. 

The Remington Rand Organization is to be thanked for, as well as 
congratulated upon, the first volume of its Instruction Manual (it is, 
however, bulky and expensive). It not only provides an excellent means 
for learning to code on the UNIVAC, but also contains enough informa- 
tion of a general nature to benefit all devotees of high-speed computa- 
tion. The text is exposed clearly enough so that a veritable tyro in the 
coding field can learn its principles solely from a perusal of this book 
without any other aid. 

The first eight chapters, which contain the essential information 
about coding on the UNIVAC, can be mastered in about two weeks by 
a beginner. Chapter 9, which gives an elementary description of the 
engineering principles involved in the operation of a high-speed com- 
puter, will prove quite interesting for many readers, although it is not 
essential for the acquisition of coding skill. This reviewer is tempted 
to make the same statement anent Chapter 10, dealing with the Flow 
Chart (since she has never found it necessary to use one), but it is sure 
to arouse the ire of the numerous advocates of this supposed aid in 
programming. 

The last Chapter serves as an Appendix and will be found most 
useful by the coder, especially as it lists all of the codes in alphabetical 
order along with complete explanations, whereas in the previous 
chapters they were given in small sets only and grouped according to 
operation. The textbook is further enhanced by the inclusion of many 
exercises in each chapter and of several excellent charts—some bound 
into the volume, others placed loose within a pocket on the inner 
back cover. 

The beginner will regret the absence of answers to the exercises with 
which he might compare the results of his first timid attempts at coding. 
He will welcome, however, the second volume dealing with advanced 
techniques in programming and coding, promised by the publisher. This 
reviewer would welcome even more a detailed operators’ manual, so 
that an advanced programmer could learn to use the machine without 
having to rely on the services of another person. 

A word of advice to a newly established laboratory planning to make 
this manual its sole textbook: One cannot expect a manufacturer to 
highlight those peculiarities of his machine which for the unwary user 
become costly and irritable pitfalls—although the UNIVAC is far less 
vulnerable on this score than are most of the other computers. It is 
imperative, therefore, to append another chapter to the manual, ex- 
pounding the danger points and means for avoiding trouble. One illustra- 
tion should suffice. The UNIVAC command for shifting to the right is 
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a negative number—a fact which must be borne firmly in mind when 
attempting to modify this command. 


I. RHODES 
NBS 


Betz, B. K., An ORDVAC Floating Binary System. Ballistic Research 
Laboratories, Aberdeen Proving Ground. 

Memorandum Report No. 733. 
NAVAL ProvinG GrounD, Dahlgren, Va. Routines for the MARK III 
Electronic Calculator. 


. MASSACHUSETTS INSTITUTE OF TECHNOLOGY. Mechanical Translation. 


Devoted to the Translation of Languages with the Aid of Machines. Vol. 1, 
No. 1, March 1954. 

The editors of this publication plan to publish it only occasionally 
and at irregular intervals in 1954. The first issue is devoted primarily 
to a bibliography of existing literature on the subject. 

Roy GOLDFINGER, New York University Compiler System. New York 
University, A.E.C. Computing Facility. Feb. 10, 1954. 


. DicitaL CoMPUTER LABORATORY, University of Illinois. JILLIAC Pro- 


gramming, A Guide to the Preparation of Problems for Solution by the 
University of Illinois Digital Computer. 

The ILLIAC is available for use by personnel of the University of 
Illinois provided users code their own problems. This manual is designed 
to teach how to code for the ILLIAC and how to make use of all its 
facilities. It includes chapters on: the arithmetic of the ILLIAC, the 
order code, subroutines and interpretive routines, the decimal order 
input, scaling, machine methods and coding tricks, checking methods, 
tape preparation, calculation of running time, preparation of a complete 
program, the cathode ray tube display, the program library, and defini- 
tion of terms. 

This, like the UNIVAC Instruction Manual and WHIRLWIND 
Comprehensive System Manual, is a good example of a text for a course 
in machine programming or a manual for self taught users of a specific 
electronic digital computer. 


J. WEGSTEIN 
NBSCL 


DiciTAL COMPUTER LABORATORY, Massachusetts Institute of Tech- 
nology. WHIRLWIND Memoranda M-2527, M-2539. 

The WHIRLWIND, operating in a supervisory mode, is used to 
simulate different easy-to-use digital computers. One of these pseudo- 
computers is called the ‘‘CS Computer’ (Comprehensive System) and 
another is called the “Summer Session Computer.”” These memoranda 
include manuals which are suitable as texts for a course in programming 
for the respective computers. These memoranda also provide a descrip- 
tion of other features of the very interesting experiments in systems of 
high speed computation being conducted at M.I.T. 
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BIBLIOGRAPHY Z 


1133. ANnon., “New digital computers developed,” Aviation Week, April 12, 

1954, p. 60-63. 

A brief account of the new Bendix Computer and ElectroData Corp.’s 
Model 203 (ElectroData Corp.—a subsidiary of Consolidated Engineer- 
ing Corp.). 

The Bendix Computer is a digital differential analyzer, operating in a 
decimal numbers system. It can be used to solve linear and non-linear differ- 
ential equations or simultaneous sets of such equations, or for integral equa- 
tions, split-boundary value problems, and individual or simultaneous sets 
of linear or nonlinear algebraic and transcendental equations. The machine 
can be programmed from punched tape or a manual control panel ; output is 
recorded on punch tape. 

The Model 203 is a high speed general purpose computer, operating 
serial fashion in binary coded decimal number system. It includes a magnetic 
drum as well as an accessory magnetic tape storage unit. Its input may be 
from standard punch-card equipment or from a specially designed photo- 
electric reader operating from previously punched tape. Output can be 
recorded on standard punch cards, tape, or be tabulated on an electric 
typewriter. 

Both computers as described by their manufacturers are in the medium 
price range. 

A. R. Cock 
NBSCL 


1134. Anon., “These machines think,’ Signal, Nov.—Dec. 1953, p. 14-16. 

This article describes the first calculator of comparable capacity to be 
produced in quantity—the IBM Data Processing Machine, the 701 which 
rents for $11,900 a month, or more, depending on storage capacity. 

Its memory devices consist of cathode ray tubes, magnetic drums and 
magnetic tapes. It is capable of performing more than 16,000 addition or 
subtraction operations a second and more than 2000 multiplication or divi- 
sion operations a second. The magnetic Drum Storage Units will store more 
than 80,000 digits and the Magnetic Tapes have a capacity of more than 
2,000,000 digits a tape. Results can be printed by means of a 150 line-a- 
minute printer at the rate of 1050 ten-digit numbers a minute or punched 
on cards at the rate of 2400 ten-digit numbers a minute or stored on mag- 
netic tapes at the rate of 1200 ten-digit numbers a second. 


R. K. ANDERSON 
NBSCL 


1135. B. V. BowbeEn, ed., Faster than thought, a symposium on digital com- 
puting machines. London, Pitman and Sons, Ltd., 1953, 414 p. Price 
35 shillings. 

Faster than thought is an interesting account of the development and 
application of electronic digital computing machines in England. In this 
symposium, Dr. BowDEN has edited the contributions of twenty-four British 
experts, including himself, to produce an entertaining and satisfying de- 
scription of the genesis of modern digital computing machines, their theory 
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and design and their multi-faceted applications in industry, commerce and 
scientific research. 

Beginning with a foreword by the Earl of HALsBury, managing director 
of the National Research Development Corporation and a preface by the 
editor, the book is divided into three main parts: Part One, The History 
and Theory of Computing Machines; Part Two, Electronic Computing 
Machines in Britain and America; Part Three, Applications of Electronic 
Computing Machines. Included also are a brief index, glossary and two 
appendices, one a reprint of pages 666-731 of TAYLor’s Scientific Memoirs, 
v. III, including editorial notes by the translator, the Countess of LOVELACE, 
the other extracts from four letters from the Lovelace Papers. 

Some indication of the balance and flavour of the book can be gleaned 
from the enumeration of the contents of its twenty-six chapters, five of 
which are included in Part One, nine in Part Two and the remaining twelve 
in Part Three. Part One includes the following: A Brief History of Computa- 
tion, The Circuit Components of Digital Computers, The Organization of a 
Typical Machine, The Construction, Performance and Maintenance of 
Digital Computers, Programming for High-speed Digital Calculating 
Machines. Part Two treats the University of Manchester Computing 
Machine, calculating machine development at Cambridge, automatic 
computation at the National Physical Laboratory, the Harwell electronic 
digital computer, the Telecommunications Research Establishment parallel 
electronic digital computer, the Imperial College Computing engine, the 
Royal Aircraft Establishment sequence-controlled calculator, calculating 
machines at the Birbeck College Computation Laboratory and six pages on 
computers in America. In Part Three, special-purpose automatic computers 
are defined and there are discussed, in separate chapters, applications of 
electronic computing machines in logical problems, crystallography, meter- 
ology, ballistics engineering, government calculations, business and com- 
merce, economics, dynamical astronomy, the playing of games—nim, 
draughts and chess. Finally, in the closing chapter there is a stimulatingly 
temperate comparison of thought and machine processes. 

With respect to Faster than thought, the reviewer is inclined to agree with 
the statement in the foreword that “It is undoubtedly the best general 
account yet written.”” It contains much material which will interest both 
specialists in the development of electronic computing machinery and 
potential users of this amazing equipment. Its account of the work of 
CHARLES BABBAGE a hundred years ago, the interest of Lady LOVELACE in 
his attempts to build large-scale automatic digital machinery and her 
pioneering work in programming, the competent discussion of computer 
applications ranging from playing chess or draughts to payroll computations 
and scientific calculations—all these are among many features of the book 
which will make it stimulating reading for the general reader and spe- 
cialist alike. 

Accomplishment in the electronic digital computer field has reached a 
magnitude which makes it impracticable to achieve completeness in a 
reasonable length exposition on the subject. It nevertheless is difficult to 
understand the omission of mention of the epoch-making large-scale Bell 
Telephone Laboratories’ Relay Computers and the parallel-mode, electro- 
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static-store National Bureau of Standards Western Automatic Computer in 
Faster than thought. Again brief mention could have been made of the versa- 
tile E. U. Conpon’s relay machine to play the game of nim demonstrated 
about two decades ago at the New York City World’s Fair, evoking con- 
siderable interest there. Finally, some indication of the trend of activity 
concerning electromagnetic and electronic digital computers on the Conti- 
nent, including the development of computers and their application, would 
not have been unwelcome. In the evaluation of the book as a whole, these 
points appear of minor consequence and Dr. BowDEN is to be congratu- 
lated upon his contribution to the literature on modern electronic digital 
computers. 


E. W. C. 


1136. D. T. Caminer, L. FAntL, G. S. D. Kine, T. R. Toompson, & W. B. 
Wricut, “The use of a high speed automatic calculator in the refinement 
stages of crystal-structure determinations,” Acta Crystallographica, v. 7, 
Part 3, 1954. 


Refinement can be said to begin when approximate atomic positions 
have been found, and the phases of the most important structure factors 
are known. Best-known methods available for refinement are 1. The Fourier 
Method, 2. The Differential Fourier Method, 3. The (Fo — F,.) Fourier 
Method and 4. The Method of Steepest Descent and Least Squares. For 
automatic calculation, the authors found method (2) most appropriate for 
the first stages of the refinement of coordinates and method (4) for the final 
stages. The Quarichi (1949) modification of (4) was used. 

The calculator used was the LEO, a machine built for the Lyons Elec- 
tronic Office. LEO is a serial, single-address, binary machine. It has a store 
of sixty-four mercury delay lines, providing a capacity of 2048 words of 17 
binary digits. It has paper-tape input, photoelectrically read. 

Programs were worked out for nicotinamide, which has the space group 
p2:/a. The computational procedure could be readily modified for other 
centrosymmetric space groups. Application to acentric space groups would 
be more complicated. 

Flow diagrams are given for the two programs of computation. An 
analysis is made of storage space required and the time needed for the 
calculations. 

E. W. C. 


1137. W. C. Carter & M. Ettts, “A comparison of order structures for 
automatic digital computers,’’ Operations Research Soc. of Amer. Jn., 
v. 2, no. 1, 1954, p. 41-58. 


A method is outlined whereby order structures of the 1-, 2-, 3-, and 
4-address types may be analyzed for application in a particular type of 
computation. The purpose is to assist in determining which will provide 
the greatest speed of operation and the greatest efficiency of storage in a 
proposed machine design. The method is applied to a payroll calculation, 
and the resultant computation times, on comparable machines with single- 
address and 3-address order structures, are presented. Some details of the 
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different types of orders assumed to be used with the two different order 
structures are given. 


R. A. Kirscu 
NBSEC 


1138. M. J. MENDELSON, “The decimal digital differential analyzer,” 
Aeronautical Engineering Review, v. 13, 1954, p. 42-54. 


This is a very clear and simple introduction to the logic of the digital 
differential analyzer. It explains the algorithm for numerical integration, 
the method of representing negative numbers and zero or negative rates, 
the organization of digital integrator registers on a magnetic drum, some 
examples of interconnections to solve simple differential equations, the rules 
governing scaling of magnitudes, and some special techniques. The latter 
include the decision integrator, the signum function, constant multipliers, 
multiplication, and division. 

R. D. ELBoURN 
NBSEC 


NEws 


The five hundred and second meeting of the American Mathematical Society was held 
at the University of Chicago on Thursday, Friday, and Saturday, April 29, 30 and May 
1, 1954. 

The program included a symposium in Applied Mathematics (sponsored by the Society 
and the Office of Ordnance Research) of which the following papers are relevant here. 


Session II W. W. LEuTERT, Ballistic Research Laboratories, Chairman. 


Computational methods, M. R. HEsTENEs, University of California, Los Angeles 

Motivations for working in numerical analysis, JoHN Topp, National Bureau of 
Standards 

Some numerical computations in ordnance problems, A. A. BENNETT, Brown University 

Western Computer Conference and Exhibit. The Western Computer Conference on 
“Trends in Computers: Automatic Control and Data Processing’’ was held at the Ambassa- 
dor Hotel in Los Angeles, California, February 11-12, 1954, under the joint sponsorship of 
the American Institute of Electrical Engineers, the Institute of Radio Engineers and the 
Association for Computing Machinery. 


Opening Ceremonies—Thursday, February 11, 10:30 AM 


Welcome Address D. H. Leumer, Professor of Mathematics, 
University of California, Berkeley 
Keynote Address W. W. McDowELL, Director of Engineering, 


International Business Machines Corp., 
Endicott, N. Y. 


Session I—Thursday, February 11, 2:00 PM 

Automatic Control—Systems Joun M. Sauzer, Hughes Research and De- 
velopment Laboratories, Culver City, 
California, Chairman 

An Experimental Digital Flight Control Mater Marco ts, J. B. Rea Company, Inc. 

System and Eric WEIss, Dynalysis, Inc. 

The DIGITAC Airborne Control System FE. E. Botites, D. W. BursBeck, W. E. 
Frapy, E. M. GrapsBe, Hughes Research 
and Development Laboratories 
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Application of Operational Digital Tech- 
niques to Industrial Control 

Digital-Analog Machine Tool Control 
System 


Experiments with a Digital Computer in 
a Simple Control System 


Session II—-Thursday, February 11, 2:00 PM 
Data Processing—Systems 


The Automatic Handling of Business 
Data 
Business Data Processing: 
A Case Study 
A. The Retailer’s Unit Control Prob- 
lem 
B. Unit Control Systems Engineering 


C. A Solution for Automatic Unit Con- 
trol 
D. The System in Operation 


Session III—Friday, February 12, 2:00 PM 
Automatic Control—Equipment 


Approaches to Design Problems in Con- 
version Equipment 


A Multi-Channel Analog-Digital Conver- 
sion System for D-C Voltages 

A High-Speed Multichannel Analog- 
Digital Converter 

A Shaft-to-Digital Encoder 


Real-Time Digital Differential Analyzer 
(“Dart”’) 


Session IV—Friday, February 12, 2:00 PM 
Data Processing—Equipment 


The IBM Magnetic Drum Calculator 
Type 650-Engineering and Design Con- 
siderations 

Remington Rand Speed Tally—Design 
Features 

Production Control with the Elecom 125 
—Electronic Business System 


BERNARD M. Gorpon, Laboratory for Elec- 
tronics, Inc. 

Harry W. MERGLER, Lewis Flight Propul- 
sion Laboratory of the National Advisory 
Committee for Aeronautics, Cleveland, 
Ohio 

T. J. Burns, J. D. Croup, J. M. Sauzer, 
Hughes Research and Development Lab- 
oratories 


RIcHARD G. CANNING, University of Cali- 
fornia at Los Angeles, Chairman 

OLIVER Wuitey, Stanford Research Insti- 
tute 


SAMUEL J. SHAFFER, The May Company, 
Los Angeles 

RaymonD Davis, Librascope, Inc., Glen- 
dale, California 

Harry D. Huskey, Institute for Numerical 
Analysis 

Myron J. MENDELSON, Computer Research 
Corporation 


RoGER Sisson, Computer Research Cor- 
poration, Hawthorne, California, Chair- 
man 

A. K. Sussk1np, Servomechanism Labora- 
tory, Massachusetts Institute of Tech- 
nology 

W. S. SHockency, Hughes Research and 
Development Laboratories 

James M. MITcHELL, J. B. Rea Company, 
Inc. 

B. M. Gorpon, M. A. Mever, R. N. 
NicoLa, Laboratory for Electronics, Inc., 
Boston, Massachusetts 

LorEN P. MEtssner, U. S. Naval Ordnance 
Laboratory, Corona, California 


Harry E. Burke, Jr., Consolidated En- 
gineering Corporation, Pasadena, Cali- 
fornia, Chairman 

E. S. HuGuHes, Jr., Endicott Laboratory 
International Business Machines Corp., 
Endicott, New York 

J. L. Hit, Engineering Research Associates, 
St. Paul, Minnesota 

NorMAN GRIESER, Underwood Corporation, 
Long Island City, New York 
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A Centralized Data Processing Center J. J. Dover, Air Force Flight Test Center, 
Edwards Air Force Base, Edwards, Cali- 
fornia 

A Merchandise Control System WiiiaM L. Martin, Telecomputing Cor- 


poration, Burbank, California 


Discussion groups Friday morning, February 12, included the following topics: 
Unit Control in Retail Operations 
Numerical Control of Petroleum and Chemical Processes 
Numerical Control of Machine Tools 
Maintenance Requirements for Business Computers 
Mathematical Methods in Management Programming 
Symposium on Automatic Programming for Digital Computers. Sponsored by the 
Navy Mathematical Computing Advisory Panel, a symposium was held in Washington, 
D. C. on Thursday and Friday, May 13 and 14, 1954. 


May 13, Morning Session: C. V. L. Smita, ONR, Washington, Chairman 
Definitions, GRacE Hoprer, Remington Rand Inc., Philadelphia 
Differentiators, HARRY KAHRIMANIAN, Naval Aviation Supply Office, Philadelphia 
Compiler Method of Automatic Programming, NorA Moser, Army Map Service, 
Washington 
Editing Generators, JOHN WaAITE, RCA Victor, Camden 
New York University Compiler System, Roy GOLDFINGER, NYU 
Application of Automatic Coding to Logical Processes, BETTY HOLBERTON, David Taylor 
Model Basin 
Discussion 
May 13, Afternoon Session: ALBERT E. SmitH, BuShips, Washington, Chairman 
The MIT Systems: Comprehensive, Summer Session, and Algebraic, C. W. ApAms and 
J. H. Lanine, Jr., MIT 
Interpretive Routines in the Illiac Library, Davin E. MULLER, Univ. of III. 
Planning Universal Semi-Automatic Coding, SAUL Gorn, Ballistics Research Lab., 
Aberdeen Proving Ground 
Present Status of the Michigan Automatic General Integrated Computation (MAGIC), 
J. H. Brown and J. W. Carr III, Univ. of Mich. 


Discussion 
May 14, Morning Session: S. B. WitLiams, Consulting Electrical Engineer, Washington, 
Chairman 
Automatic Programming on the Burroughs Computer, HUBERT M. LIvinGsTON, Burroughs 
Corporation 


IBM 701 Speed Coding System, Joun W. Backus and HaRLAN HERRICK, International 
Business Machines Inc., New York 

Programming for the IBM 701 with Repetitively Used Functions, ALLEN KELLER and 
R. A. BuTTERWortTH, General Electric Company, Lynn, Mass. 

Discussion 

Summary and Forecast, GRACE Hopper, Remington Rand Inc. 

National Applied Mathematics Laboratories. The National Bureau of Standards Insti- 
tute for Numerical Analysis, operated on the Los Angeles Campus of the University of 
California since 1947, ceased to exist on June 30, 1954. Most of its activities will be continued 
under the sponsorship of the Office of Naval Research and the Office of Ordnance Research 
under the administration of the Department of Mathematics of the University of California 
at Los Angeles of which M. R. HEsTENEs is chairman; C. B. TompKins will be the Director 
of the new unit which will be called Numerical Analysis Research, [NAR]. The physical 
facilities of INA, including the SWAC and the library, will be loaned by the National 
Bureau of Standards to the University of California. 
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A new section, Numerical Analysis, has been established July 1, 1954 by the National 
Bureau of Standards in Washington; at the same time, the Machine Development Labora- 
tory has been replaced by a section Mathematical Physics. 

The new organization of the National Applied Mathematics Laboratories will be as 
follows: 

F. L. Aut, Acting Chief, E. W. Cannon, Assistant Chief 
. Numerical Analysis, J. Topp, Chief 
. Computation Laboratory, M. ABRamowi7z, Acting Chief 
. Statistical Engineering Laboratory, C. EisEnnarRt, Chief 
. Mathematical Physics, E. W. Cannon, Chief 


Pwr 
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BIBLIOGRAPHY Z 


1139. R. R. BENNETT, “Analog computing applied to noise studies,” I. R. E., 
Proc., v. 41, 1953, p. 1509-1513. 


A problem of great importance in modern communication engineering 
is the analysis of noise in linear and nonlinear systems. For a linear system 
the weighting function (impulsive response, Green’s function) completely 
characterizes the system and is a necessary tool in spectrum shaping 
analysis. The author discusses the analog simulation of the weighting func- 
tion using the concept of the adjoint system (cf. H. H. LANNtnG, Jr. & R. 
H. Battin, MTAC, v. 7, 1953, p. 125). For nonlinear systems the 
author briefly describes how the amplitude distribution may be shaped 
using a function generator. Finally, statistical methods—using the concept 
of confidence levels—are discussed for determining the accuracy, say of the 
standard deviation, of a number of identical runs on the same physical 
system. 

K. S. MILLER 
New York University 
New York, N. Y. 


1140. M. B. Coy ze, “An airflow analogy for the solution of transient heat 
conduction problem,” British Jn. Applied Physics, v. 2, 1951, p. 12-17. 


The airflow analogy parallels closely the familiar electrical resistance- 
capacity analogy for transient heat conduction phenomena. Quantity of heat 
and temperature in the heat conduction system are represented in the airflow 
system by mass of air and pressure respectively. The continuous thermal 
system is considered replaced by a lumped system which is simulated by a 
network of air reservoirs, each containing a mass of air which is a function 
of the air pressure within it, inter-connected by capillary tubes through 
which the mass airflow is proportional to the pressure difference between 
the ends. The air reservoirs are in the form of U-tubes with a confining fluid, 
all the U-tubes having one leg in common. With this arrangement, the 
capacity (taken in the electrical sense to mean the mass of air discharged 
for unit fall in pressure) of a reservoir may be made a function of pressure 
by constructing it with variable cross-sectional area. Therefore, a problem 
involving variable thermal properties may be set up with appropriately 
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constructed reservoirs and a solution obtained with no further manipulation 
required during the test. 

This appears to be an advantage of the air-flow analogy relative to the 
electrical analogy since in the latter variable thermal properties are simulated 
by varying circuit elements in the course of a test, either manually or auto- 
matically. However, the electrical analogy appears to possess far greater 
flexibility and inherent experimental accuracy. 

J. H. WEINER 
Columbia University 
New York, N. Y. 


1141. LEONARD J. CraicG, “The magnetic amplifier as an analog computer 
component,” I. R. E., Proc., v. 41, 1953, p. 1477-1482. 


It is shown that the load voltage of a parallel-connected magnetic 
amplifier is a gate whose width is directly proportional to the control cur- 
rent, when a square wave source voltage is applied. Tests indicate that this 
relationship holds to within 0.1 percent over a wide range of source voltage 
(except for scale factor), load resistance, and temperature, provided the 
magnetic material used has reasonably ideal characteristics. 

This phenomenon is utilized in several circuits: a linear pulse width 
modulator, a linear pulse-position generator and a square-wave phase shifter. 
The phase shifter was applied to a DC multiplier which achieved product 
accuracies of 3 percent, without the use of high precision components. 


Joun B. WALSH 
Department of Electrical Engineering 
Columbia University 
New York, N. Y. 


1142. W. A. McCoot, ‘An AM-FM electronic analog multiplier,” I. R. E., 
Proc., v. 41, 1953, p. 1470-1477. 


A multiplier is described wherein an AM signal whose index of modula- 
tion is proportional to an input voltage is applied to an attenuator whose 
gain is proportional to frequency deviation (an FM discriminator). Fre- 
quency deviation is made proportional to a second input, so the discriminator 
output, with the carrier subtracted, is the product of the inputs. Feedback 
techniques stabilize the AM and FM modulation indices. Although basically 
a two-quadrant multiplier, it is arranged, by a conventional method, for 
four-quadrant operation. 

A total of 24 tube envelopes is used. Performance is limited primarily by 
stability of discriminator parameters, and is independent of tube character- 
istics. A long time multiplication error of less than 1 percent of maximum 
output is reported. Usable bandwidth (flat amplitude response and less 
than 15 degree phase shift for highest frequency component) is approxi- 
mately 0.1 percent of carrier frequency (1 percent was reported, apparently 
erroneously). 

Joun B. WALSH 
Department of Electrical Engineering 
Columbia University 
New York, N. Y. 
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1143. James C. Ross, “A calculator for aiding matrix calculations,” 
Faraday Soc. Trans., v. 50, 1954, p. 8-12. 


Twelve Curta Calculators are mounted on a base and have their working 
mechanisms geared together and centrally controlled. The operations (1) 
shift carriage, (2) clear and (3) enter multiplier, activate the twelve calcu- 
lators simultaneously. Complete constructional details are given. The 
capacities are, of course, those of the Curta—6D multiplier, 8D multiplicand 
and 11D product. Provision for positive and negative multiplier and multi- 
plicand is possible. 

The device is well suited to problems involving many simultaneous 
multiplications and addition. As an example, matrix multiplication is de- 
scribed. The writer claims the equipment is suitable for use with inversion 
schemes like those of Crout, but no further explanation is given. For twelfth 
order matrix by vector multiplication with 4D elements, a time of 17 
minutes is cited. 

A. OPLER 
The Dow Chemical Company 
Pittsburg, Calif. 


1144. O. J. M. Situ, “Economic analogs,’”’ I. R. E., Proc., v. 41, 1953, 
p. 1514-1519. 


The possibility of using electronic analog computers for problems in 
economics is discussed. Several circuits are given for specific simple problems 
and economic terms are interpreted in terms of electrical networks; for 
example, if the supply (outcome) influences the demand (input) then there 
is feedback in the system. The author points out that the concepts of closed 
loop system, Nyquist diagrams and vector diagrams are now being used by 
economists to interpret and analyze a variety of economic problems. A list 
of forty references is included. 

K. S. MILLER 
New York University 
New York, N. Y. 


1145. G. W. Swenson & T. J. Hicomns, “A direct-current network analyzer 
for solving wave-equation boundary-value problems,” Jn. Appl. Phys., 
v. 23, 1952, p. 126-127. 


The authors discuss the solution of the partial differential equation 
Ve + (w/c?e = 0 


on the function g(x, y) by means of a resistance network. This network is 
analogous to the usual network for the Laplacian, except that the term 
(w/c)? requires the introduction of a “negative resistance” involving an 
amplifier between each node and ground. This “negative resistance’ must 
be readjusted manually by an iterative method in the course of the problem. 
Both characteristic value problems and forcing function problems are 
considered. 

F. J. M. 
Columbia University 
New York, N. Y. 
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1146. G. H. VinEyarp, ‘Simulation of trajectories of charge particles in 
magnetic fields,’ Jn. Appl. Phys., v. 23, 1952, p. 35-39. 


The author sets up an equivalence between the motion of a charged 
particle in an electric and magnetic field and the motion of a sphere rolling 
on a rotating surface, subject to gravity and viewed from an independently 
rotating coordinate system. A variety of interesting special cases corre- 
sponding to the betatron and magnetron are discussed. In addition an experi- 
mental technique involving a rotating glass surface and multi-exposure 
photography for obtaining the orbits is described. A preliminary setup was 
made with immediately available components and the author considered 
the results satisfactory. 

F. J. M. 
Columbia University 
New York, N. Y. 


NOTES 


163.—ANOTHER UNRECORDED Pitiscus ITEM. In my article on Bar- 
THOLOMAUs PiTiscus (1561-1613), to whom our word Trigonometry is due 
(1595) (MTAC, v. 3, p. 390-397, 1949), I referred to one of his tables, 
“unlisted in any of the ordinary bibliographic or historical sources”’ (p. 396), 
entitled 


SINUUM/TANGENTIUM /et/SECANTIUM/Canon Manualis/Acco- 

modatus ad trigo/nometriam/BARTOLOMAEI/Pitisci/Griinbergensis 

Silesij./[Decoration ]/HEIDELBERGAE,/Typis Iohan Lincelloti, Acad 

Typo./Imprensis Iondé Rosdée./MDCXIII./Signatures A”-H®, I‘ 

[200 p.]. 

It contained 5D tables for all six of the trigonometric functions in the 
quadrant, at interval 1’. Opposite pages record the values for 30 minutes of 
each function; but both pages are headed “‘sinus,”’ “‘tangens,’’ ‘‘secans.”’ 
The angles on the right-hand pages being complementary to those on the 
left, the values on the right-hand pages are really those of cos, cot, csc for 
the angles on the left-hand page. Thus 180 pages, A? verso—H® recto, are 
occupied with the tables. ‘‘Explicatio numerorum hujus Canonis” occupies 
pages A! verso to A? recto. A! recto is the title page. H® verso—-H™ recto are 
devoted to text explanations; H” verso and I! verso blank; I' recto comment 
on the following 230 Errata (sin and cos 46; tan and cot 86; sec and csc 96; 
2 others) : I? recto—I* recto. 

The little pages are of size 7.3 X 13.2 cm. A film copy of this table for 
the Library of Brown University was made from a copy of the original at 
the University of Illinois. 

In the summer of 1953 the bookshop Old Authors Farm, R.R. 1, Harris- 
burg, Ontario, Canada, offered an extraordinary collection of mint copies of 
old books—duplicates from the Vatican Library. Among these was a second 
edition of the Pitiscus volume described above. The displays of the title 
pages, down through the word HEIDELBERGAE, except for a new decora- 
tion, are identical; then follows in the new volume: Typis Joan. Georg. 
Geyder. Acad. Typ. Imprensis Jondé Rosdé./MDCXX./A"-H®. This 
volume, acquired by Mr. ALBERT E. Lownes of Providence, R. I., further 
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enhancing the value of his remarkable personal library in the field of history 
of science, was kindly placed at my disposal for preparing this article. 

It will be observed that there are eight fewer pages in the new edition; 
this was due to elimination of the 230 errata noted in the first edition, and 
correction of the corresponding places in the earlier pages. I checked the 46 
noted errors for sines and cosines and found that while one of those was 
incorrectly listed, not only were all the others corrected, but also some errors 
not previously noted. In spot checking of the other functions no error was 
found in the new edition. Otherwise pages A*-H” in both editions appear to 
be textually identical. 

The little volume is bound in old vellum, seemingly as fresh as on the 
day the volume was published 334 years ago—seven years after the death 
of Pitiscus. 

Apart from this copy, the only other one known is in the Vatican Library. 

R. C. ARCHIBALD 
Brown University 
Providence, R. I. 


164.—PERTURBATION OF LAGRANGIAN COEFFICIENTS. Suppose it is de- 
sired to interpolate a function g(x) polynomially to a certain value when the 
independent variable x is a function, f(t), of a parameter, ¢, where the values 
g are known at equally spaced values of ¢. Under certain circumstances we 
find it advantageous to carry out a first order perturbation of the Lagrangian 
coefficients for equally spaced values. Suppose f(¢) is a function with a non- 
vanishing derivative and with a continuous third derivative. Also, sup- 
pose we want to interpolate g(x) to the value x = f(a), given f(t) for 
t=a+t (i — b)h, t = 0, ---, m where b is a positive number less than m 
and h is a positive number. Then, from Taylor’s series 


(1) f(a+ @ — b)h) = f(a) + @ — d)hf'(a) + 3 — bYhF" (a) + R: 


where R; is the appropriate remainder. For the purpose of interpolation 
f(a) has no effect on the coefficients. Since the values (¢ — b)hf'(a), 4 = 0, 
+++, m are equally spaced, the Lagrangian coefficients for interpolating to 
x are given by 

(e — x) 


(2) L; - 1 


j = 0,1, ---, m. 
fae 7 


Now, considering the x’s as variables we find the logarithmic derivative of 


L; to be 
dL; dx — dx; dx; — dx; 
3 dL; _ ye — dey 
* L; tem yy — 
Then letting 
= — d)hf'(a), dx; = 3 — bP HF"(a) 


x = 0, dx = 0 
dL; _ hf'(a) i) — 
ee. 





= 2f’ (a) [m(b — j)] j = 0,---,m 








184 NOTES 


Hence the perturbed values of the Lagrangian coefficients are given by 
b-—i hf"’ (a) : | 
(4) Ly + ats = 1 (F=3) [1+ FO mo —a]. 


The sum of the perturbed coefficients in (4) equals 1 since 


E[o-om(t=)]-» 


While an error estimate might be calculated, in any application it is simple 
to calculate actual coefficients to compare with the first order approximations. 

For example, let f(t) = #, ¢: = a+ (¢ — 1.5)h toi = 0, 1, 2, 3. Then 
f(a) = 6a, f’(a) = 3a? and for interpolation to x = a? from 4 











1 Oh 
Lo + dLo = - 36(1 +3) 
9 3h 
n+dh= 7 (1+3) 


Lz + dLz 


= 


$f, . #4) 
6 2a 


1 Oh 
ke +dle = = a(! -3). 


h 
Clearly these coefficients should be used only if r is small. Formulas of this 


type when appropriate have been used at the Los Alamos Scientific Labora- 
tory in a calculation on the IBM 701 to save appreciable storage space 
compared to that required for direct calculation of the correct coefficients. 
Of course, f’’(a)/f’(a@) may not be simple to calculate in many applications. 

In case it is inconvenient to calculate f’(a) and f(a) we may interpolate 
to x = f(a) using the following adjusted coefficients 


(6) Ly+dL; =] (22> )fi + a m(b +b -a| 





j = 0, ---,m, 


where ¢ is a convenient value of t, usually taken to be one of the equally 
spaced values near a and p = (a — c)/h. When a = ¢, (5) reduces to (4). 
The sum of the adjusted coefficients in (5) is 1. 

PRESTON C. HAMMER 
Univ. of Wisconsin 
Madison, Wis. 

Harwoop G. KoLsky 
Los Alamos Scientific Lab. 
Los Alamos, N. M. 


165.—AN INVESTIGATION INTO THE REAL Roots OF CERTAIN POLY- 
NOMIALS. In connection with work on the application of statistical theory 
to the determination of the number and location of roots of high order poly- 
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nomial equations, it was required to investigate the real roots of the 
equations: 


(1) Cte’ t+ xe8t+ eo ett Ht HPtxt+1=0 


where all possible combinations of sign are taken. 
It can be shown simply, that all roots of these equations lie in the range: 


4 < |x| <2 


and the computational problem therefore reduces to that of evaluating each 
of the 256 different equations obtained by changing the signs of (1), over 
the range $ < x < 2, and at a sufficiently fine interval to ensure that no 
roots are missed. 

Preliminary calculations indicated that the functions are quite well 
behaved, and that tabulation at intervals of ;#; in x would be adequate, and 
it was thus necessary to evaluate each equation at 25 points. 

The polynomials were written in the form: 


x(+++x(x(x +1) +1)---) 41 


and were evaluated by successive addition of the coefficients and multiplica- 
tion by x. 
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186 NOTES 


The calculations were performed on the HEC 1 and the time for calcu- 
lation and printing of each polynomial was 50 seconds. For convenience of 
reference, the sign combination of each polynomial was printed as well. 

An analysis of the results of this tabulation shows that of the 256 equa- 
tions considered : 

58 have no real roots 
190 have two real roots 
8 have four real roots. 


The distribution of the roots by numerical value is shown in Fig. 1. 
The detailed results of this tabulation are available for reference at the 
Computation Laboratory, Birkbeck College. 

The author wishes to thank Mr. R. L. MicHAELson of the British Tabu- 
lating Machine Company for providing computing facilities for this 
calculation. 

K. H. V. Bootu 
Computation Laboratory 
Birkbeck College 
21, Torrington Square 
London, W.C. 1. 


166.—NUMERICAL STUDY OF SIGNATURE RANK OF CuBIC CYCLOTOMIC 
Units. In the search for algebraic fields (particularly beyond the quadratic) 
where unique factorization fails, little work seems feasible on high speed 
computers for the reason that the straightforward evaluation of forms must 
be controlled by group-theoretic data-processing.! A rare opportunity, 
however, is afforded by the cubic field R(a) generated by the Gaussian sum 


p—1 
a= > exp 2rix*/p 
z=0 
where » = 6m + 1 isa prime. Here the so-called cyclotomic unit? is given by 
the simple formula (for three conjugates 0 < i < 2): 


n—1 
@, = — [I sin (Sug'**4/p) osc (Jng'41/p) 
t=0 


where g is an odd primitive root modulo p (or an even one augmented by #). 
Now while the decimal accuracy of such a formula would be lost very rapidly 
the 0; have the property that the possible non-unique factorization in R(a) 
can be of a certain frequent type (namely even class number) only when all 
three 0; are positive.? This condition is easy to determine by machine. 

We accordingly form three tallies Ao, Ai, Az for each p as follows: We 
increment the A; tally (0 <7 < 2) when the least positive residues of 
g**t#+1 and g** modulo 2,0 < t < m — 1, lie in different halves of the interval 
(0, 2p) subdivided at p. Such a procedure tallies the possible negative sign 
of sin (}rg**+*+1/p) csc (42g**/p); hence the three 0; are seen to be all 
positive exactly when the Ao, Ai, and Az have (final) values that are odd, 
even, and odd respectively. 

Since there are many primitive roots g for each p we find that the A; 
depend on the g chosen, except for the fact that the occurrence of three 
positive 6; must be independent? of g. Assuming that the powers of g lay 
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in half-intervals “‘at random’’ we should expect the A; to be each approxi- 
mately p/12, while assuming “random parities,”’ (which is more precarious), 
we should expect the 6; to be all positive approximately } of the time; 
(since 699,02 = + 1, or, equivalently, A» is always odd). 

The computation was performed on the MIDAC starting Jan. 27, 1954, 
for the first 207 values of p (<3000). The values* of » and g in the form of 
binary-coded decimals were placed as alternate words on high-speed photo- 
electric input tape. They were read in by pairs (taking approximately .15 
sec. for each pair). The conversion to binary form was internal. For the main 
induction loop the starting point was each g**, 0 < ¢ < m — 1, from which 
were formed g*‘t*+! for 0 < i < 2 (always modulo 2p). The three tally 
increment decisions (for A;) were made as described above (taking .04 sec. 
for each circuit of the loop). In addition, every time the MIDAC chose not 
to increment an A; tally it incremented a common D tally. After the loop 
had been traversed n times two checks on accuracy were made. First the 
sum of the three A; and the D tallies was verified to be (p— 1)/2 (the 
number of increment decisions), and secondly the last power, g-»/? (re- 
duced modulo 2), was verified to be 2 — 1. The MIDAC was made to 
examine the three final A; tallies for the occurrence of alternate parities 
and indicate such occurrences or non-occurrences in the form of a — 0 or 0 
respectively in a temporary storage space. The output for each 0 consisted of 


pb g Ao Ai Az — 0 (or 0) 


where the first five items were converted to decimal internally and shifted 
for a short-word (4 character) print-out while the sixth item was the signed 
zero indicating presence or absence of alternate parities. (Print-out time 
for each p was approximately 5 seconds.) 

The total running time was about one hour of which 20 minutes was 
input-output. 

The 21 values of (<3000) for which the A; have alternate parity (or 
for which the 6; are all positive) are reproduced below with j their order in 
the list of primes (of the form 6” + 1). 


j Pp j p j p 

18 163 71 853 143 1879 
27 277 77 937 145 1951 
33 349 81 1009 156 2131 
37 397 107 1399 169 2311 
48 547 129 1699 191 2689 
52 607 135 1777 199 2797 
60 709 137 1789 200 2803 


It will be seen that the frequency of these p is close to #5, which is not in 
agreement with ‘‘randomness”’ of parities (as postulated above). The magni- 
tudes of the A; did come out, however, to agree rather well with the random- 
ness assumption made earlier, and they are omitted here. A spot check‘ of 
pb = 163, 277, 349, 2803 by hand calculations revealed that unique factoriza- 
tion does fail in the field R(a), and that the class number is 4 for each of 
these cases. 


188 NOTES 


The author is indebted to Dr. J. W. Carr III and the staff of the Willow 
Run Research Center for kindly making the MIDAC available for this study 
(which is part of Army Ordnance Project DA-20-018-ORD-12332). 


HARVEY COHN 
Department of Mathematics 
Wayne University 
Detroit, Michigan 
1 Compare D. H. ot, Guide to tables in the theory of numbers. National Research 
Council, 1941, p. 75-77, TAussky, Some ns, problems in algebraic number 
theory. National Bureau FS Standards Report (to appea: 


? For theoretical background consult H. Hasse, y= OO Bestimmung von Grundein- 
heit. Berlin, 1950, p. 70. 
* These were taken from a table of minimum positive g for p < 3000 in I. M. Vino- 


GRADOV, Osnovy Teorii Chisel [Fundamentals of the Theory of Numbers]. Moscow, 1940, p. 
110. 


* The case p = 163 was discovered through another procedure by E. Artin, according 
to a private communication. 


167.—CULLEN NuMBERsS. These are numbers of the form 2" + 1 and 
are remarkable in that they seem to be composite for m > 1, although there 
is no @ priori reason for this. CUNNINGHAM & WOODALL! made a study of 
these numbers and found them all composite with a small factor for 
1 <n < 141. No factor of 141-2'4* + 1 is known. I have completely 
factored the following cases left incomplete by Cunningham. The case 
n = 46 is due to R. A. Lrénarp of Lyons. 


n2" + 1 n n2" + 1 


47 -6031230671 42 23-43-83-2250270487 

37 - 32502455213 43 3-5-163-2633-58752797 
3-5-339016085231 45 11-47-2437-1256655529 

3?- 20879 - 55586743 46 5-31-47-139297-3189821 

41 - 3433 - 152326961 47 7-11-43-3593-556021079 
41-131611-8150491 48 7-379-997-5107973329 

13-43-1291 - 124932557 66 5*-13-67-107- 131-8353 - 382030403 


N. G. W. H. BEEGER 
Nicolaas Witsenkade 10 
Amsterdam 


1A. J. C. CunnincHam & H. J. Woopatt, “Factorisation of Q = (24 gq) and 
(q:2¢ + 1),” Messenger Math., v. 47, 1917, p. 1-38. 


CORRIGENDA 


V. 6, p. 225, 1. 11, for monomial read elementary. 
V. 7, p. 34, |. 6, for 6 read 1. 
V. 7, p. 175, 1. 17, for 9 read 8. 








